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ABSTRACT 

A  theory  and  resultant  calculations  are  presented,  predicting 
the  noise  generated  by  a  lifting  rotor  in  forward  flight.  With 
suitable  adjustments  in  parameters,  this  theory  reproduced 
(by  numerical  integration)  the  closed  form  results  of  Garrick 
and  Watkins  in  their  investigation  of  noise  produced  by  a  pro¬ 
peller.  The  theory  presented  here  accounts  for  (1)  asymmetry 
of  inflow,  lift,  ana  drag;  (2)  non-linear  section  characteris¬ 
tics,  compressibility,  and  reverse  flow;  (3)  first  harmonic 
and  steady  rigid  blade  flapping  and  pitching;  (4)  rotor  shaft 
tilt;  and  (5)  chordwise  distribution  of  lift  and  drag  forces. 
Sound  pressure  levels  were  evaluated  for  up  to  20  harmonics  at 
a  general  field  point  translating  with  the  rotor  hub,  for  the 
H-34  and  HU-1  helicopters.  A  comparison  with  experiment  is 
also  presented. 
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SUMMARY 


The  lifting  rotor  is  considered  as  a  swept  surface,  area 
segments  of  which  are  subjected  to  oscillating  pressures, 
expressed  as  a  Fourier  series  in  time,  as  was  done  for  the 
propeller  in  axial  flow  by  Garrick  and  Watkins.5  In  this 
case,  of  course,  the  theory  is  extended  tc  include  in-plane 
components  of  forward  speed  and  azimuthal  asymmetry.  The 
sound  pressure  at  any  field  point  is  then  found  by  a  straight¬ 
forward  numerical  integration  which  uses  as  inputs  (1)  lift 
per  unit  span  as  a  function  of  radius  and  azimuth,  (2)  for¬ 
ward  speed,  (3)  rate  of  descent  or  climb,  (4)  RPM,  (5)  rotor 
geometry,  including  non-linear  compressible  airfoil  section 
characteristics,  and  (6)  first  harmonic  and  steady  blade  root 
angles  of  flapping  and  pitching.  This  approach  has  the  ad¬ 
vantage  of  calculating  either  near  or  far  field  noise,  includ¬ 
ing  all  effects  other  than  "thickness  pressure"5  and  that  due 
to  the  "white"  components  attributable  to  random  turbulence 
in  the  boundary  layer  or  the  wake.  Impulsive  changes  in  load 
around  the  azimuth,  due  either  to  compressibility  or  to  vor¬ 
tex  interactions  are  therefore  easily  investigated.  Radiation 
patterns  are  investigated  for  the  H-34  and  HU-1A  helicopters, 
both  around  the  azimuth  and  in  the  shaft  direction,  for 
several  forward  speeds,  and  the  sound  spectra  compared  with 
measurements  for  selected  cases. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


The  theory  for  calculating  propeller  rotational  noise  has 
been  extended  to  account  for  (1)  asymmetry  of  inflow,  lift, 
and  drag;  (2)  non-linear  section  characteristics,  compressi¬ 
bility,  and  reverse  flow;  (3)  first  harmonic  and  steady  rigid 
blade  flapping  and  pitching;  (4)  rotor  shaft  tilt;  and  (5) 
chordwise  distribution  of  lift  and  drag  forces,  such  as  are 
found  in  the  operation  of  lifting  rotor-propellers  in  for¬ 
ward  flight.  A  numerical  integration  is  used  to  evaluate  the 
sound  pressure  levels  for  up  to  20  harmonics  at  a  general 
field  point  translating  with  the  rotor  hub. 

By  comparison  with  a  closed  form  solution,  it  has  been  shown 
that  the  numerical  integration  is  quite  accurate  for  axial 
flow  when  azimuthal  increments  of  1°  (or  smaller)  are  used. 

In  addition,  the  pronounced  reduction  in  first  harmonic  sound 
level  indicated  by  the  closed  form  solution  at  a  flight  Mach 
number  of  two-tenths,  two  diameters  to  the  side  and  about 
three-fourths  of  a  diameter  in  front  of  the  propellers5,  is 
shown  to  be  a  spurious  consequence  of  a  "far  field"  assumption. 
Similar  nodes  at  other  flight  speeds  and  field  points  would 
probably  also  be  severely  modified  by  an  "exact"  numerical 
integration  using  the  subject  theory. 

Application  to  two  helicopters  currently  in  widespread  use 
showed  that  the  first  harmonic  rotational  noise  has  by  far 
the  highest  sound  pressure  levels  in  hover,  but  because  of 
the  low  frequencies  is  still  marginal  with  respect  to  the 
threshold  of  hearing.  The  second  harmonic  frequency  falls 
just  above  the  lowest  for  which  there  is  experimental  data, 
and  there  is  good  correlation  with  measured  sound  pressure 
levels  for  this  component.  Successively  higher  harmonics, 
however,  are  predicted  to  fall  off  rapidly,  whereas  the 
measured  values  'do  not.  Contaminating  influences  in  the 
test  setup  (which  were  not  performed  for  the  present  pur¬ 
pose)  include  the  presence  of  the  tail  rotor,  reciprocating 
engine  exhaust  noise,  "sounding-board"  effect  of  fuselage 
and  other  surfa'ces,  and  relatively  large  uncertainty  as  to 
the  elevation  of  the  hub  with  respect  to  the  microphone.  De¬ 
ficiencies  in  the  theory  include  lack  of  "thickness"  rota¬ 
tional  noise,  "vortex"  noise  due  to  boundary  layer  turbulence, 
wake  vortex  noise,  effect  of  ground  plane  and  downwash  field 
refraction  on  the  radiation  pattern,  and  frequency  shifts 
due  to  time-varying  velocities  in  the  downwash  field. 

The  principal  effect  of  forward  speed  is  to  raise  the  higher 
harmonic  levels.  This  increase  may  occur  abruptly  or  gradu¬ 
ally,  depending  on  the  harmonic  content  of  the  aerodynamic 
section  loadings,  as  would  be  expected.  This  increase  occurs 
much  more  gradually  with  the  four-bladed  H-34  than  with  the 
two-bladed  HU-1.  No  attempt  at  generalizations  should  be 
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made  without  more  extenaive  inveatigation  than  waa  poaaible 
in  thia  study. 

Octave  band  radiation  patterns  in  the  azimuthal  sense  were 
close  to  symmetric ,  particularly  for  the  HU-1  helicopter,  one 
diameter  beneath  the  hub  over  a  speed  range  from  hover  to 
115  knots  and  at  radii  of  two  and  ten  diameters.  Some  small 
polar  dissymmetry  appears  in  the  azimuthal  radiation  patterns 
of  the  H-34  helicopter,  for  the  lower  octave  band,  at  the 
higher  speed,  and  at  ten  radial  diameters. 

Variations  in  octave  band  sound  pressure  level  were  noted  with 
change  in  vertical  distance  beneath  the  rotor.  The  variations 
are  similar  to  those  for  the  first  harmohic  of  a  propeller  in 
axial  flow.  For  the  rotor,  a  peak  occurs  between  two  and  five 
diameters  below,  depending  upon  radial  distance  from  the  hub. 

The  reduction  of  sound  pressure  level  with  radius  in  the  first 
three  octave  bands,  at  an  azimuth  ^O*270#,  is  roughly  as  one 
would  expect.  However,  if  the  onset  of  spherical  spreading  is 
chosen  as  the  criterion  for  determining  "far  field",  this  may 
be  more  like  one  hundred  diameters  from  the  hub  than  three,  as 
is  often  used.  Moreover,  a  purely  mathematical  definition  of 
far  field  may  lead  to  spurious  theoretical  predictions,  as 
shown  by  the  propeller  example  mentioned  above. 

Increasing  the  number  of  blades  appears  to  be  the  most  effec¬ 
tive  means  of  reducing  rotor  rotational  noise.  An  example  in 
which  the  number  of  blades  was  increased  from  two  to  three  for 
the  HU-1  resulted  in  a  significant  decrease  in  sound  pressure 
level.  The  effects  of  compressibility,  RPM  and  impulsive  down- 
wash  were  not  significant.  Increasing  blade  chord  resulted  in 
a  small  but  fairly  general  decrease  in  sound  pressure  level. 

The  following  appear  to  be  promising  directions  for  future 
works 

1.  Perform  sound  level  measurements  using  existing  rotors 
under  conditions  specifically  designed  to  isolate  rotor 
noise  from  all  other  sources.  (Some  of  this  work  has  un¬ 
doubtedly  already  been  done  in  research  carried  out  at 
Southampton  and  perhaps  United  Aircraft,  but  data  have  not 
been  generally  available.)  Correlate  the  levels  so  ob¬ 
tained  with  the  results  of  the  present  theory  as  applied 
to  the  new  test  configurations. 

2.  Using  analytical  techniques  such  as  proposed  in  Ref.  10, 
and  perhaps  measurements  on  rotating  cylinders  and  in  the 
boundary  layer  on  full-scale  rotors,  attempt  to  predict 
vortex  noise  for  correlation  with  Refs.  21  and  22,  and  ul¬ 
timately  arrive  at  a  method  for  calculating  rotor  vortex 
noise  to  be  added  to  the  rotational  noise  of  the  present 
study. 
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3.  Extend  the  present  digital  computer  program  to  allow 
prediction  of  the  rotational  sound  fields  of  tandem 
rotors  and  rotor r  in  ground  effect. 

4.  Attempt  to  predict  the  effect  on  rotor  noise  of  both 
ordered  and  random  velocities  in  the  wake  beneath  the 
rotor . 
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INTRODUCTION 


Understanding  the  lifting  rotor-propeller  as  a  noise  source 
has  recently  become  a  technical  objective  of  major  importance. 
There  are  two  factors:  (1)  military  use  of  helicopters  in 
tactical  situations  imposes  the  requirements  for  low  external 
noise  levels  with  a  new  authority1;  and  (2)  the  introduction 
of  turbine  powerplants  in  the  present  generation  of  rotary¬ 
wing  and  V/STOL  aircraft  leaves  the  rotor-propeller  as  the 
dominant  source  of  low  frequency  noise.  These  components, 
of  course,  are  precisely  those  which  are  attenuated  least 
with  distance.  Thus,  the  aircraft's  rotor-propeller  noise 
characteristics  will,  in  combination  with  its  approach  speed, 
determine  how  much  time  an  enemy  will  have  to  prepare  for 
the  aircraft's  arrival. 

A  second  important  consideration  is  the  difficulty  of  attain¬ 
ing  a  satisfactory  sound  treatment  to  reduce  rotor  noise  as 
experienced  inside  the  aircraft.  It  is  generally  accepted 
(Ref.  2,  Figures  9  and  10,  for  example)  that  such  reductions 
must  be  achieved  at  the  source,  for  frequencies  below  about 
200  cps,  since  the  weight  penalty  for  significant  attenuation 
is  usually  prohibitive. 

From  both  these  viewpoints,  as  well  as  consideration  of  the 
obvious  influences  of  rotor  noise  on  helicopter  ASW  operations, 
commercial  applications,  and  the  relatively  unexplored  opera¬ 
tion  of  VTOL  propellers  in  the  transition  regime,  it  is  im¬ 
portant  that  greater  insight  be  achieved  regarding  the  mecha¬ 
nism  of  rotor-propeller  noise  generation. 

Investigations  in  the  field  of  helicopter  rotor  noise  have  the 
benefit  of  years  of  work  dealing  with  airplane  propellers. 3  * 4  • 5 
Differences  between  measured  rotor  noise  levels  and  those 
predicted,  however,  using  the  empirical  relations  for  propel¬ 
lers6  indicate  that  the  applicable  range  of  parameters  is  too 
limited,  even  when  used  for  hover  or  vertical  flight.  The 
profound  differences  which  come  about  in  forward  flight  are, 
of  course,  clear.  It  is  the  latter  to  which  this  study 
specifically  addressed  itself. 

The  theory  of  propellers  is  extended  here  to  allow  prediction 
of  the  rotational  noise  generated  by  lifting  rotors  in  forward 
flight.  It  will  be  assumed  that  the  time-varying  pressures 
experienced  at  each  radial  station  along  a  blade  are  known  for 
a  given  helicopter  configuration  and  operating  condition.  From 
this,  the  sound  radiation  amplitude  and  distribution  around  the 
rotor'  azimuth  will  be  predicted  in  the  fixed,  that  is,  non¬ 
rotating  system.  With  these  results,  it  will  be  possible  to 
treat  the  rotor  like  any  other  noise  source  on  the  aircraft  in 
predicting  what  is  experienced  in  the  aircraft  or  by  am  ob¬ 
server  on  the  ground. 
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The  effects  of  the  approximately  "white"  noise  which  exists 
in  the  turbulent  boundary  layer  sheathing  the  blades,  and  is 
modulated  for  the  observer  by  the  blade  passage  frequencies7, 
will  be  ignored  in  this  study.  Also,  pressure  fluctuations 
experienced  as  a  result  of  the  passage  of  the  complex  wake 
of  distributed  shed  and  trailed  vorticity  with  resoect  to  a 
non-rotating  observer  translating  with  the  rotor8#9  will  not 
be  considered.  Although  some  recent  analyses10  promise  sim¬ 
plification,  the  approach  of,  and  in  fact  the  extension  sug¬ 
gested  in,  Ref.  5  is  followed  rather  directly.  Effects 
specifically  accounted  for  in  the  sections  to  follow  include: 

1.  Variations  in  velocity  normal  to  the  leading  edge,  angle 

of  attack,  lift  and  drag  with  radius  and  around  the  azimuth. 

2.  Non-linear  and  compressible  airfoil  section  characteris¬ 
tics. 

3.  The  influence  of  non-uniform  induced  flow  on  angle  of 
attack,  and  (in  those  instances  where  steady,  two- 
dimensional  section  data  indicate  that  the  section  should 
be  stalled,  yet  it  is  not)  on  velocity  normal  to  the 
leading  edge. 

4.  Reversed  flow. 

5.  Coning,  flapping,  and  cyclic  pitch  variations. 

The  effects  of  rotor-fuselage  or  rotor-rotor  interference 
will  not  be  specifically  examined  except  as  they  may  exist 
in  the  experimental  or  theoretical  spanwise  loadings  used  as 
input  for  the  analysis.  It  is  felt,  however,  that  the  theo¬ 
retical  investigation  reported  here  is  a  necessary  first  step 
in  increasing  understanding  of  lifting  rotor-propeller  noise. 

The  underlying  philosophy  in  this  approach  stems  from  the 
impulsive  nature  of  the  airloads  on  a  lifting  rotor  in  for¬ 
ward  flight.  Rabbott  and  Churchill11  first  showed  the  very 
rapid  fluctuations  which  could  be  expected  with  azimuthal  po¬ 
sition;  Miller8,  Piziali  and  DuWaldt9  and  others  tied  the 
existence  of  these  very  high  harmonics  theoretically  to  the 
proximity  of  vortex  elements  trailed  from  the  passage  of  pre¬ 
ceding  blades  or  the  same  blade  in  preceding  revolutions. 

Aside  from  the  fact  that  approximations  for  other  quantities  - 
such  as  amplitude  and  phase  radius5  -  would  probably  have  to 
be  made  for  a  closed  form  integration  over  the  propeller- 
rotor  disc,  it  was  felt  that  analytical  representation  of  the 
harmonics  of  airloads  associated  with  non-uniform  downwash 
would  be  likely  to  lose  components  of  interest.  Since  numeri¬ 
cal  integration  would  obviate  the  need  for  both  kinds  of 
approximations,  this  approach  was  taken.  The  problems  of 
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obtaining  satisfactory  numerical  accuracy  and  limiting 
digital  computer  machine  time  to  reasonable  amounts  were, 
of  course,  recognized  from  the  outset,  and  will  be  discussed 
in  the  Results  section  of  this  report. 


ANALYSIS 


This  formulation  follows  the  general  procedure  developed  in 
Ref.  5  and  preserves  much  of  the  notation,  the  exceptions 
occurring  where  it  conflicts  with  helicopter  terminology. 

COORDINATE  SYSTEM 

The  coordinate  system  is  Cartesian,  rectangular,  right- 
handed  and  oriented  so  that  (1)  the  z-axis  is  vertical,  (2) 
the  positive  x-direction  points  in  the  direction  of  flight, 
and  (3)  the  origin  is  taken  at  the  center  of  the  rotor  hub 
(see  Figure  1) ,  so  that  the  system  translates  with  the  flight 
velocity,  V. 

Points  on  the  surface  area  swept  out  by  the  rotor  blades  will 
be  called  "source  points"  and  identified  as  points  x,  y,  z. 

It  will  be  convenient  to  use  polar  coordinates  on  the  area 
swept  by  the  rotor  blades.  The  azimuth  angle,  ijs,  will  be  de¬ 
fined  in  a  plane  perpendicular  to  the  shaft  and  taken  to  be 
zero  on  the  negative  x-axis,  as  shown  in  Figure  2. 

A  general  point  outside  the  swept  area  but  translating  with 
the  rotor  will  be  called  a  "field  point"  and  denoted  by  xQ, 
yQ,  zQ.  Defining  field  point  azimuth  relative  to  the  pilot, 
the  azimuth  angle  of  a  field  point,  ty  ,  is  such  that  ^  =0 
on  the  positive  x-axis.  Since  an  integration  will  eventually 
be  performed  over  the  swept  area,  the  "source  points"  will  be 
specified  using  a  double  subscript  notation  ir  j.  Here  i 
identifies  the  radius  position  and  j  the  azimuthal  position. 

Equations  for  determining  x,  y,  z  coordinates  of  the  source 
points  are  derived  geometrically  by  examining  the  position  of 
a  particular  blade  of  the  rotpr  at  any  time,  t.  It  is  assumed 
that  the  blade  is  inflexible  but  has  a  "flapping  hinge",  off¬ 
set  a  distance,  e,  from  the  center  of  rotation.  This  allows 
a  blade  to  rotate  out  of  the  plane  of  rotation  through  an  angle 
6^,  which  is  assumed  to  be  small.  A  point  on  the  blade  is  de¬ 
fined  as  being  a  distance  r  from  the  hub  when  0^*0.  No  other 

motion  of  the  blades  with  respect  to  the  hub,  other  than  pitch 
changes  (denoted  by  e^j) ,  will  be  accounted  for.  The  rotor 

shaft,  however,  is  pictured  as  tilted  through  an  angle  a  in 

the  longitudinal  plane.  Blade  point  coordinates  are  then  given 
(see,  for  example,  Ref.  12),  when  small  angle  assumptions  are 
made  for  sin  6  and  cos  0,  by 
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Fig. 2.  COORDINATE  SYSTEM. 
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THE  ACOUSTIC  EQUATION 
As  in  Ref.  5,  we  define 

p<xo'Vzo't)AVVVV*o,ei“tl 


-  WWV1"80*1 


as  the  sinusoidally  varying  pressure  at  a  field  point  (x  , 
yo,  *Q)  with  the  frequency  w-mBfl.  ° 

Here  fl  is  the  rotational  frequency  of  the  rotor,  and  R ,  of 
course,  signifies  taking  the  real  part  of  a  complex  number. 
Steady  conditions  will  be  assumed)  i.e. ,  what  happens  in  one 
revolution  happens  in  every  other,  and  in  fact,  what  happens 
to  one  blade  at  a  particular  azimuth  is  repeated  on  every 
other  blade  when  it  is  at  that  azimuth.  It  follows  that  all 
phenomena  will  be  periodic  with  the  fundamental  period 
2 

T  ■  rj  thus  m  is  the  harmonic  number.  If  we  can  determine 

the  harmonic  components  of  the  forces  acting  on  the  fluid  at 

a  source  point  in  the  directions  of  the  coordinate  axes,  namely 

,  and  F  ,  then  we  can  calculate  the  sound  pres- 
m^  m  m 

*ij  Yij  zij 

sure  at  a  field  point  from  the  expression 


R«  lpm*xo'yo'zo* 


-iko 


n  rl— {p  _L_+f  i  +p  }e  imBflt. 

RelTTFm  Sx^ra  JyFm  Tz^) § e  1 
x  o  y  Jo  z  o 


(2) 


which  is  given,  in  slightly  different  form,  as  Eq.  (6)  in 
Ref.  5.  Here, 
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mBfl 


S  -  /(xQ-x) 2+y*  I(y0-y)2  +  (z0-z)2] 
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M(xq-x)+S 

y2 


C  ■  velocity  of  sound 

y  -  /I^mT  (2a) 


REPRESENTATION  OF  FLUID  FORCES  IN  TERMS  OF  FORCES  ON  THE 
BLADES 


The  forces  on  the  fluid  are  those  which  cause  the  rotor  to  ex¬ 
perience,  as  reactions,  lift  and  drag.  An  area  element  of  the 
swept  surface  may  be  thought  of  as  having  such  forces  on  it  as 
arise  from  calculating  a  force  per  unit  of  projected  blade  area 
and  then  multiplying  by  the  differential  of  swept  area.  This 
step  may  be  obvious,  but  if  it  is  kept  in  mind,  it  will  em¬ 
phasize  that  we  are  dealing  with  "effective  areas"  and  "effec¬ 
tive  pressures."  Losing  sight  of  this  might  make  it  appear 
that  we  are  calculating  fluid  pressures  in  the  x-direction 
different  from  those  in  the  y-direction,  or,  as  another 
example,  that  we  are  thinking  of  viscous  drag  forces  as  arising 
from  fluid  pressures  just  as  normal  forces  do;  in  fact,  neither 
of  these  statements  is  true. 


Let  us  examine  the  force  per  dnit  blade  area,  irrespective  of 
its  source  or  direction,  which  is  experienced  by  an  element  of 
the  swept  surface  at  some  azimuthal  and  radial  station.  It 
will  have  a  time  history  such  that  as  each  rotor  blade  sweeps 
through,  pulses  are  experienced.  The  duration  of  these  pul¬ 
ses  is  Tij*  and  they  recur  after  a  period  T,  as  shown  in 

Figure  3.  Note  that  if  the  blade  chord  length  is  and  the 


blade  geometric  pitch  angle  0, .,  then  the  projection  of  the 
chord  length  on  the  swept  surface  is  C^cose^  and  ■ 


C^cose 

”"FPr 


To  avoid  the  complications  which  otherwise  arise 


i 


it  will  be  assumed  from  this  point  that 


Cicos6ij*Cicos0i4C^  and 


where  0^  is  the  mean  pitch  angle  at  a  given  radius.  Since 


X  ~6 

eij  “  e0+9icos^j)  *92gin  j)  "R-e 


where  g  is  the  twist  built  into  the  blade  in  degrees  per 


11 


I  ? m-*r*w+* 


gjfirnTr  —  -t  . 
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fraction  of  blade  radius ,  approximating  0^  by  0^  is  tanta¬ 
mount  to  neglecting  61  and  02/  which  are  rarely  higher  than 
about  10° ,  and  the  total  angle  0,.,  will  almost  always  be  less 
than  25°.  13 

The  shape  of  the  pulses  will,  of  course,  depend  on  the  distri¬ 
bution  of  force  over  the  chord.  Now  suppose  we  limit  ourselves 
to  operating  regimes  where  the  rotor  blade  is  un stalled  and  re¬ 
verse  flow  can  be  ignored.  If  we  handle  lift,  incompressible 
drag  and  wave  drag  separately,  then  we  can  assume  with  little 
error  that  only  the  magnitude  of  the  pulses  in  Figure  3  will 
vary  with  azimuth,  and  that  their  shapes  will  be  constant, 
except  for  the  variation  of  t:  with  radius.  In  Ref.  5,  the 
effect  of  chordwise  force  distribution  is  discussed  at  some 
length  as  having  appreciable  effects  only  on  the  higher  har¬ 
monics.  On  the  other  hand,  the  rotor  "bang"  or  "slap"  noise 
reported  in  Refs.  7  and  13  is  a  relatively  high  frequency 
phenomenon,  and  is  such  an  important  aspect  of  rotor  noise 
that  chordwise  distribution  differences  have  been  purposely 
retained  in  this  analysis.  Putting  aside  for  the  moment  the 
question  of  the  magnitude  of  the  pressure  pulses,  their  shapes 
can  be  expressed  as  a  real  Fourier  series; 


a_  +  cosmBflt+b 

O  .  *»  _  m .  1 

i  ro 


■lmi 


m. 


sinmBflt) 


Note  that  the  shape  function,  f. ,  can  be  thought  of  as  a  func¬ 
tion  either  of  chord  length  projection  on  the  swept  surface,  C! 
or  of  time,  t.  This  function  lacks  the  index  j,  identifying  1 
azimuth,  since  it  is  thoughtof  as  "normalized,"  and,  also  for 
the  moment,  as  beginning  at  t*0,  which  coincides  with  i|>j*0.  For 

convenience  the  coefficients  a_  ,  a_  ,  and  b.  are  normalized 

so  that  f.  *1.0.  Thus,  to  account  for  the  variation  in  lift 
max 

and  drag  around  the  azimuth,  we  must  establish  normalization 
factors  for  these  Fourier  coefficients. 

Since  we  have  specified  that  the  shape  represented  by  the  nor¬ 
malized  coefficients  will  always  have  a  maximum  ordinate  equal 
to  unity,  then  the  aerodynamic  lift  acting  on  a  span  length 
Ar^  is 

Lij(t)  “  CL,  .^(Uij,2ciAriA*ij{fidciAri 
13  chord 
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where 


iT^  ■  normalization  factor  for  lift 

U, .  ■  resultant  relative  wind  velocity  normal  to  blade 

3  leading  edge 

p  ■  air  mass  density 


lift  coefficient 


Similarly,  for  incompressible  profile  drag 

»V  <«  *  CD  ^<V2VriA"ij|VCi“i 


chord 


and  for  compressible  (wave)  drag 


Dw.  j  (t>  *CD„  ^(uij)2ci4riAKSij|iiaq1ri 


chord 


Since  sound  is  an  oscillatory  phenomenon,  we  may  drop  the 

steady  term  in  the  Fourier  series.  Furthermore,  since  we  have 

expressed  the  amplitudes  of  the  force  components  F  ,  F  , 

mx  mv 

ij  yij 

and  Fm  as  the  complex  multiplier  of  the  time  function 

zij 

^imBnt 

e 

we  will  prefer  to  rewrite  the  real  Fourier  series  in  the  form 

00  00 

7,  *  cosmBftt+b_  sinmBftt) -Vc  cos(mBnt-f  .) 

1  m*lmi  mi  m«lmi 


Therefore,  the  normalized  coefficient  c  e  mi  is  associated 

mi 

with  the  blade  lift,  by  the  equation  directly  above;  then  we 
can  define  two  more  sets  of  similar  quantities  associated  with 
(1)  incompressible  drag,  say  d  and  r  .,  and  (2)  compressi- 
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bility  drag,  namely  and 

Analytical  expressions  for  such  Fourier  constants  are  de¬ 
rived  by  the  usual  means  (see,  for  example.  Ref.  14) ,  for  the 
three  chordwise  shapes  shown  in  Figures  4(a)  through  (c) .  The 
selection  of  these  shapes  as  approximations  for  lift,  incom¬ 
pressible  drag  and  wave  drag  was  based  on  a  perusal  of  two- 
dimensional  section  data  below  stall,  and  computational  con¬ 
venience. 

The  corresponding  formulae  are  as  follows: 
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Fig. 4.  IDEALIZED  FORCE  DISTRIBUTION  VERSUS  CHORD. 
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Note  that  the  lift,  incompressible  drag  and  wave  drag  co¬ 
efficients  are  all  functions  of  (1)  local  "free  stream"  Mach 
number,  ,  perpendicular  to  the  blade  leading  edge  and  (2) 

aerodynamic  angle  of  attack  a^. 

To  allow  for  the  fact  that  these  pulses  do  not  always  start 
at  t*0,  but  at  time  t+tj/j/n,  we  need  merely  multiply  by 


-imBft  (♦  ./n)  -imBij/ , 
e  j  "  e  3 


where 

♦j  4  ( j“l) A* 

A 1 1>  •  increment  in  azimuth  angle  from  one  station  on  the 
swept  area  to  another. 


j  “1,  2,  3  •  •  •  n 


n  ■  total  number  of  azimuthal  sections  into  which  the 
swept  area  is  divided 

From  these  definitions  it  follows  that  the  m^*1  harmonic  com¬ 
ponent  of  the  forces  experienced  by  the  fluid  due  to  blade  lift 
and  drag  may  be  expressed  as  minus  the  following: 


17 


R  [L 
e  iQj 


imBnt, 
e  ] 


y»cL  |un]2ci 

t  ij  i  .-if_J.-imB*J.imB(ltt_  . 

■  R  { — r  ■  c  e  mie  rje  Jr.A^Ar. 

•  IrTdc! "i  1  1 


J  i  i 
chord 


,  eimB0tj 
e 


D  d 

vij  mi  -ir. 


Re(y>ciIuijl2|-r=i^re‘  r,ni  + 


jfidC- 

chord 


'D  e 
wii  ni.-iK 

wr 

chord 


jein’Bnt)ri4i(i4ri 


Before  these  sinusoidally  varying  forces  can  be  evaluated  or 
expressed  in  terms  of  components  along  the  coordinate  axes, 
it  is  necessary  to  establish  (1)  the  direction  and  magnitude 
of  the  resultant  velocities  U.  .  (since  lift  and  drag  forces 
are  defined  as  perpendicular  and  parallel  to  them,  respec¬ 
tively)  ,  and  (2)  the  blade  orientation. 

RESULTANT  VELOCITY  AT  A  BLADE  SECTION 

Rotor  aerodynamics  derivations  (see,  for  example,  Ref.  12) 
usually  start  by  defining  two  velocity  components  in  a  plane 
normal  to  the  blade  leading  edge.  These  are  Up  and  U^/  com¬ 
ponents  parallel  and  perpendicular  to  the  shaft  axis,  respec¬ 
tively.  The  convention  for  Up  is  positive  up;  UT,  of  course, 

is  always  taken  positive  when  directed  toward  the  airfoil 
leading  edge. 

Accounting  for  components  of  forward  speed,  rotational  speed, 
an  average  induced  velocity  component  v  (assumed  positive  down) 
and  an  effective  cgmponent  normal  to  the  chord  as  a  result  of 
flapping  velocity  6,  yields 
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U_  ■  Vcosa  sim|/+fir 

X  S 

U  *  V(8ina  -Bcosa  cos^i)  -v-B  (r-e) 

P  S3 


This  neglects  the  radial  component  of  velocity,  assumes  that 
each  strip  of  airfoil  perpendicular  to  the  leading  edge  acts 
as  an  airfoil  in  two-dimensional  flow  with  the  same  resultant 
Velocity  normal  to  the  leading  edge  and  the  same  angle  of 
attack,  and  makes  the  usual  small  angle  assumptions.  The  ac¬ 
tual  induced  velocity  is  a  quantity,  which,  in  general,  varies 
both  over  the  blade  radius  and  azimuth,  and  its  calculation 
has  been  the  subject  of  extensive  research.  Its  effect  on 
angle  of  attack  may  be  quite  significant,  especially  where  U 

might  otherwise  be  small.  It  is  usually  assumed,  however,  that 
induced  velocities  in  the  plane  of  the  rotor  are  small  compared 
to  U^,  except  where  that  velocity  is  small,  in  which  case  the 

resulting  forces  are  unimportant.  For  this  reason,  the  re¬ 
sulting  velocity  will  be  calculated  using  an  induced  veloc¬ 
ity  v  assumed  to  be  uniform  over  the  disc,  and  calculated  from 
the  momentum  theory  expression  (Ref.  15) . 


v 


_ thrust _ _ 

2nR2 p [ (Vs in a  -v)  2+ (Vcosa  )2]^^ 
s  s 


The  approach  used  to  account  for  the  actual  variations  in  in¬ 
duced  velocity  over  the  disc,  both  normal  and  in  the  plane  of 
the  rotor,  will  be  explained  in  the  following  section. 

DETERMINATION  OF  BLADE  ELEMENT  ANGLE  OF  ATTACK,  LIFT  AND  DRAG 
- v - — — - '•  ■ 

As  shown  in  Figure  5,  the  blade  angle  of  attack,  ^ ,  is  a 
function  of  the  velocity  components  and  UT  and  the  local 
geometric  blade  pitch  angle,  6^. 

In  this  study,  it  will  be  assumed  that  the  lift,  or  more  ex¬ 
actly,  the  loading  per  unit  span  normal  to  the  blade  chord, 
the  "aerodynamic  section  loading,"  is  known  from  experiments 
or  previous  calculations.  Actually,  measurements  reported  in 
Ref.  16  and  17  are  to  be  used  as  inputs  to  this  analysis. 

This  eliminates  the  need  for  an  accurate  knowledge  of  v  as  a 
function  of  r  and  <ji,  which  is  acknowledged  as  a  necessity  if 
anything  but  the  lowest  harmonics  oi  blade  airloads  is  needed. 

The  procedure  is  as  follows.  Normal  force  coefficients  CN(M,0) 

will  be  tabulated  for  the  airfoil  section  of  interest,  as  will 
the  aerodynamic  section  loading  (ASL)  for  a  given  flight  con- 
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LIFT 


DRAG 


Fig. 5.  BLADE  SECTION  QUANTITIES. 
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dition  as  a  function  of  r  and  i| >.  For  each  element  of  swept 
area  corresponding  to  a  source  point,  the  resultant  velocity 
will  be  calculated  from 


*  /(Vcosa  sin^  .  +fir . )  *+(V[sina  -B  .cosa  cosi|» .  )-v-B  .  [r-e] )  2 

Sjl  SjSj  J 

ui-l 

the  Mach  number  from  M^j  ■ 

and  the  normal  force  coefficient  from 

(ASL) 


\ ^ t 

:Ni j  “  l/2p(Uij)2ciAri  * 


(5) 


The  table  of  airfoil  coefficients  then  will  be  entered  with  the 
known  values  of  CM  and  M..,  and  the  corresponding  value  of 

Nij  13 

aij  read  out,  interpolating  when  necessary.  With  and  aij 

thus  known,  similar  airfoil  coefficient  tables  for  C.  and 

Lij 


'D. 


as  functions  of  these  same  variables  will  be  entered  and 


ij 


the  lift  and  drag  coefficients  read.  A  simultaneous  " look-up' 
for  Cn  at  the  corresponding  but  with  will  be  made 


'ij 


‘ij 


to  determine  CQ 
uv 


This  subtracted  from  Cn  will  yield 

°ij 


CD  .  This  information  is  all  that  is  needed  in  addition  to 

Wij 

Eg.  (3)  and  Eq.  (5)  to  completely  determine  L  and  D  as 

mij  raij 

written  in  Eg.  (4)  . 


It  is  important  to  note  that  the  value  of  CN  given  by  Eq.  (5) 

may  exceed  any  that  can  be  found  in  the  tables  of  two-dimen¬ 
sional  coefficients,  as  pointed  out  in  Ref.  18.  This  can  be 
explained  in  several  ways:  (1)  the  in-plane  components  induced 
by  the  vortex  wake  may  be  signif icantly  increasing  the  resul¬ 
tant  velocity  ^ ,  (2)  the  effects  of  radial  flow  may  be  suf¬ 
ficiently  energizing  the  boundary  layer  to  delay  stall,  (3) 
unsteady  effects  may  be  sufficient  to  delay  stall.  In  view  of 
the  paucity  of  airfoil  section  data  either  in  the  presence  of 
spanwise  flow  or  under  unsteady  conditions,  the  approach  taken 
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■-..-  u- 

a»pf 


here  waa  to  have  the  machine  program  evaluate  the  minimum 
increment  required  in  U4  .  to  make  CM  fall  within  the  steady, 

13  Nij 

two-dimensional  data.  For  that  radius  and  azimuthal  station, 
the  adjusted  value  of  will  thereafter  be  used. 

Certain  details  of  the  "table  look-up"  operation  involving 
Cg,  CL,  CD,  and  a  are  worth  noting.  The  present  program  is 

set  up  for  symmetric  airfoils;  all  values,  therefore,  are  tab¬ 
ulated  as  positive.  The  extension  to  cambered  sections  would 
not  be  difficult.  In  any  case  the  sign  of  the  quantities  UT 

and  Up  do  not  appear  explicitly,  so  care  must  be  taken  to  dif¬ 
ferentiate  between  regions  of  positive  and  negative  angle  of 
attack,  and  between  normal  and  reversed  flow.  The  signs  of  UT 

and  (ASL)  are  tested  and  the  procedures  listed  in  Table  I  are 
followed. 

Considerations  leading  to  Eq.  (6)  in  the  section  to  follow 
will  show  that  these  special  procedures  represent  the  aero¬ 
dynamic  forces  as  vectors  in  the  appropriate  direction.  This 
is  a  straightforward  transformation  of  coordinates  operation; 
again,  using  the  approximations  cosB3l  and  sinB-B,  it  follows 
that 

F  ■  [  (Lsin+-Dcos+)  cosa_sin^+(Lcos$+Dsin$)  (Bcosa  cosi|i-sina  )  ] 

X  8  8  8 

Fy  -  [  (Lcos$+Dsin$)  Bsinij/- (Lsin<fr-Dcos$)  cos<j>] 

F  ■  [  (Lsin+-Dcos+)  sini^sina  +(Lcos$+Dsint)  (Bsina  costjj+cosa  )  ] 

Z  8  8  8 

(6) 

where,  for  convenience,  the  subscripts  m,  i  and  j  are  omit.-.ed 
and  considered  to  be  understood. 

APPLICATION  OF  THE  ACOUSTIC  EQUATION 

Now  to  use  Eq.  (2)  to  calculate  the  mth  harmonic  of  sound  pres¬ 
sure  at  a  field  point  (x  .  y  .  z  )  due  to  the  nr“  harmonic 

oscillatory  pressure  at  a  source  point  identified  by  the  in¬ 
dices  (i,j),  we  need  only  (a)  to  note  that 


22 


a  «“iJcrn0^  -iko . ,  ik  M  ik  .  (x  -x) 

3/6  m  id*  _  -e  m  id  r  m  ,  /  m,  i  *  o  , 

ay  '  <?.  .  1  fi  |,'T7+fi.  '  s.  .  J 


3xo  Sij 


’ij 


'  Si j  Sij 


a  «"ik«°<n  “ik  a ,  ,  2 

.3_(i _  m  l -e  .  u  ^  flk  t  y  1  (v  -v) 

a„  v  c  )  (S,.)*  11  m  S.  1  Ky^  y) 


Tyo  Sij 


'ij 


’H 


(2o“2) 


(7) 


(b)  to  substitute  L  and  D  from  Eq.  (4)  for  L  and  D,  re- 

mij  mij 

spectively,  in  Eq.  (6) ,  and  (c)  to  substitute  the  resulting 
force  components  and  the  Eq.  (7)  expressions  into  Eq.  (2) . 
The  resulting  equation  can  be  broken  into  real  and  imaginary 
parts  and,  after  algebraic  manipulation,  cast  in  the  follow¬ 
ing  form: 

1  Ki-t  _ 

WWij,,  "  TVTs^^ij^/^ij-^^ij1  - 


ij 


‘ij 


S^M 


ij 


B«.ij(riJSin*ij+JijC08*ij)1> 


1  _  _  _ 

VWVijj  ■  17  Ts^T5-(Aij1Dmijain»ij+c»ijC08*ijI  * 


S .  .M 


Fi»ij<Tij8in<iJ+JijC08*iJ)1> 


(8) 
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where 

Kij  "  ^2(yo-yij)co8^j-[(xo-xij)co80l8  + 

Y2  (zQ-zij)  sinas]8ini|/j } 

®ij  "  nx0-xij)jj+Y2[(y0-yij)ajsin^j  + 


(zQ-zij  )$  jCOSi(ijSinas+co8as)  ] } 
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m 


ij 


D„  d 
Vij  mi 


cos(k  o,.  +  r  +mBi|».)  + 
fidCi  "‘I  *!  1 


I 


D  e 
wii 

■m — - cos  (k_o  ,  ,+k  +mBiK) 

fidCi  m  ij  mi  j 


T«  ■  (B^jCosaacos(|;j-8inag) 


K±a  *  {g—“F  +vT®m  ^ 
ij  Sij  mi;J  yT  mij 


Eij  “ 


and 


K^j  *  l/2pC^  (U^j )  ^r^Aijiir^ 


Sij  *  / <xo'xi  j  ’  ‘+^  1  1  V^i  j  >  ‘  +  ( zo'zi  j  > 7 1 


•«  -  Mixo-y+sa . 
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The  integration  proceeds  by  summing  all  the  real  parts  of  the 
oscillatory  pressure  of  harmonic  m,  contributed  by  all  the 
combinations  of  i  and  j;  this  also  applies  to  the  imaginary 
parts  of  the  sound  pressure  components.  If  Pm  and  Pm  are 

R  I 

in  psi,  then  the  sound  pressure  level  in  decibels  is 


>2+(??V,  )2‘ 


dbm  '  10  lo^10" 


ij 


ij 


0.0002 


27 


RESULTS 


The  dangers  inherent  in  a  numerical  integration  of  an  ex¬ 
pression  of  the  sort  given  in  Eg.  (8)  lie  in  the  alternating 
nat\ire_of  the  trigonometric  functions  in  the  expressions  for 
E_  ,  F_,  ,  ,  and  H  .  The  arguments  of  these  functions 

“ij  mij  mij  mij 

vary  around  the  swept  area  of  the  rotor-propeller  as  the  sum 

(mB4/ .+k  o .  . )  .  The  quantity  k  will  be  recognized  as  the  in- 
j  m  13  m  th 

verse  of  the  wave  length  of  the  m  component  of  sound,  and 

thus  km0^j  is  the  number  of  waves  of  this  component  between 

the  source  point  at  ij  and  the  field  point.  This  fact  is  em¬ 
bodied  in  the  definition  of  0^  as  the  "phase  radius."  Clear¬ 
ly,  as  (mB^j+kmo^j)  becomes  a  larger  and  larger  number,  more 

and  more  significant  figures  must  be  carried  to  obtain  a  mean¬ 
ingful  answer,  since  only  that  part  left  after  subtracting  2tt 
times  the  largest  integer  leaving  a  positive  quantity  is  per¬ 
tinent.  The  term  mB^,  of  course,  may  be  thoughtof  as  "good" 

to  an  infinite  number  of  places,  and  since  k  =mBfl/C  this  can 

m 

also  be  thought  of  as  infinitely  "good",  any  error  being  dis¬ 
missed  as  making  exact  some  other  case  with  a  slightly  differ¬ 
ent  rotational  speed,  n.  The  phase  radius  0^,  therefore,  is 

critical.  The  digital  program  is  arranged,  accordingly,  to 
iterate  on  the  value  of  Si j ,  using  a  recurrence  scheme  where 


sijk 


So+iSk 


where 

tsiiki*-s0M*sk]* 

k+1  2 (S  +ASV) 

O  K 

The  computer  continues  to  iterate  until 


ASijJc+1  <  10”5feet. 

-4 

This,  for  example,  is  within  about  1/3x10  wave  lengths  for 
the  20th  harmonic  of  a  six-bladed  rotor  turning  at  300  RPM. 

As  a  test  both  of  the  program  itself  and  of  the  accuracy  of 
the  numerical  integration,  a  calculation  was  performed  for  the 
sound  pressure  level  at  the  fundamental  frequency  generated  by 
a  two-bladed,  10-foot-diameter  propeller,  producing  1600  pounds 
of  thrust  at  1582  RPM  and  a  forward  speed  corresponding  to 
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M« 0 . 2 ,  and  drawing  2680  foot-pounds  of  torque. 

Results  for  this  case  are  given  in  Figure  5  of  Ref.  5  as  a 
function  of  xQ  when  yQ*20  feet,  zQ»0.  The  present  program  is, 

of  course,  capable  of  the  same  calculation,  merely  by  setting 
as«90°,  eij*0o*  Bj"0,  the  remaining  input  appropriately. 

If  only  this  is  done,  however,  a  one-to-one  correspondence 
would  not  exist  between  the  Garrick  and  Watkins  results  and 
those  obtained  here.  The  following  approximations,  made  in 
the  earlier  work,  must  also  be  introduced  in  the  subject  anal¬ 
ysis,  at  least  for  this  one  check  case: 

ri  “  'effective  *  °-8R  }  *f^°tive 


Bp/apc^u^lc^ 


co8<p.  .+Cn  sin$.  .  ]Ar.  ■  thrust 
13  uij  1 

*  1600  pounds 

cos<|> .  . -CT  sin$.  .]Ar.  *  torque 
13  "i  j  13  1 

»  2680  foot-pounds 


uni¬ 

form 

lift 

and 

drag 


mBC 

“Jr 


^mi  *  Kmi  *  rmi 


0 


approximate  square  wave 


y2y  y 

S  2  SQ — s  — •  -  in  exponentials 

“  S  in  denominator 

of  fractions 

(g-) 2  neglected  compared  to  g 


1 

?  " far  field"  approximation 

J 


The  case  described  above  was  run  both  with  and  without  the 
"far  field"  and  square  wave  approximations,  using  an  interval 
of  1°.  The  results  are  compared  with  those  of  Ref.  5  in  Figure 
6.  It  is  interesting  to  note  that  (1)  the  numerical  integra¬ 
tion  duplicates  the  closed  form  analysis,  (2)  the  sharp  reduc¬ 
tion  in  sound  level  predicted  roughly  one  diameter  ahead  of  the 
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Sound  Level ,  db 


DIAMETER  •  10ft 
M-0.2  THRUST  -1,600  lb. 

y*  -  20ft  TORQUE  -  2,680  f  t.-lb 

Zo-0  B-2 


0  NUMERICAL  INTEGRATION, 
FORM  OTHERWISE  THAT 

OF  REF.  5 
♦  AS  ABOVE,  BUT 
CORRECTING  FOR  S 
APPROXIMATION 
B  AS  ABOVE,  BUT 
CORRECTING  FOR  SQUARE 
WAVE  FOURIER 

APPROXIMATION 


Normalized  Distance,  diameters 


Fig.  6.  FIRST  HARMONIC  PROPELLER 
NOISE  COMPARISON. 
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propeller  appears  to  be  spurious,  and  a  result  of  the  Nfar 
field"  approximation,  and  (3)  for  this  harmonic,  the  square 
wave  approximation  leads  to  no  significant  errors. 

COMPARISON  WITH  HELICOPTER  NOISE  MEASUREMENTS 

Although  this  correlation  with  earlier  computations  is  en¬ 
couraging,  the  crucial  test  for  any  prediction  method  is  a 
comparison  with  experimental  results.  Simultaneous  aerody¬ 
namic  section  loadings  and  sound  pressure  measurements  for  an 
isolated  rotor  were  not  available;  however,  such  measurements 
are  contained  individually  in  Refs.  16  and  19,  respectively, 
for  the  complete  H-34  helicopter  in  a  hover.  The  gross  weights, 
moreover,  were  roughly  the  same  in  the  two  tests.  According¬ 
ly  the  sound  pressure  levels  were  calculated  for  the  H-34 
helicopter  in  hover  for  220  main  rotor  RPM  at  a  radial  dis¬ 
tance  200  feet  from  the  rotor  hub,  and  at  a  field  point  azi¬ 
muth  angle  of  150°.  The  radial  and  azimuthal  coordinates  of 
the  field  point  correspond  closely  to  "Position  23"  in  Ref.  19, 
and  these  horizontal  distances  could  be  well  controlled  in  the 
referenced  tests.  Vertical  distances  were  subject  to  a 
greater  percentage  of  uncertainty,  since  the  pilot  was  attempt¬ 
ing  to  hover  the  aircraft  with  the  wheels  5  feet  off  the  ground, 
and  the  microphones  were  hand-held.  It  will  be  shown  later 
that  from  hub  height  to  two  rotor  diameters  beneath  the  rotor, 
sound  levels  in  an  octave  band  can  vary  as  much  as  25  decibels. 
Accordingly,  the  zq  distance  was  chosen  as  -10  feet  and  -16 

feet  to  indicate  the  sensitivity  of  this  parameter.  These 
dimensions  were  chosen  as,  roughly,  engine  height  and  wheel 
height,  respectively. 

The  results  of  the  calculation  are  shown  in  Figure  7  as  ver¬ 
tical  bars.  To  approximate  measurements  made  over  a  reflective 
surface,  all  the  levels  indicated  as  "calculated"  have  actu¬ 
ally  had  3  decibels  added  to  them. 

The  comparison  with  the  measured  values  is  notable  more  in 
the  lack  of  correlation  than  in  the  agreement  obtained.  The 
first  harmonic  frequency  is  below  the  response  range  of  the 
instrumentation  used,  so  no  check  can  be  made  for  this  com¬ 
ponent.  The  second  harmonic  falls  very  close  to  a  cluster  of 
the  measured  peaks,  but  for  the  third  and  higher  harmonics 
calculated  levels  are  substantially  low. 

The  source  of  discrepancy  can  be  sought  in  both  the  theory  and 
the  test  results.  Consider  first  the  latter.  The  tail  rotor 
and  the  reciprocating  engine  frequencies  are  indicated  on  the 
same  figure;  note  that  the  tail  rotor  fundamental  is  exactly 
6  times  that  of  the  main  rotor.  With  a  view  toward  more  re¬ 
fined  test  data  as  well  as  to  substantiate  the  data  presented 
in  Ref.  19,  a  second  narrow-band  analysis,  not  presented  in 
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Main  Rotor  Harmonic  Number,  m 

Fig.7.  COMPARATIVE  SOUND  LEVELS,  CALCULATED  versus  TEST. 


Ref.  19,  was  obtained.  This  second  sound  spectrum  is  for  the 
same  aircraft  under  the  same  operating  conditions,  but  with 
the  field  point  azimuth  shifted  to  i|>ob300°;  at  this  field 

point,  maximum  masking  of  the  tail  rotor  and  engine  exhaust 
noise  would  be  expected.  Some  idea  as  to  the  contribution  of 
sources  other  than  the  main  rotor  is  given  by  comparing  the 
test  curves  for  t|/o«300°  and  i|>o-150°  in  Figure  7.  Since  the 

masking  certainly  does  not  eliminate  the  effect  of  the  unwanted 

noise  sources,  "corrected"  test  data  would  be  still  lower  than 

that  for  «300°. 
o 

It  is  reasonable  to  question  the  number  of  azimuthal  stations 
needed  to  properly  represent  harmonics  as  high  as  the  twen¬ 
tieth.  This  case  was,  therefore,  run  with  various  numbers  of 
azimuthal  increments.  These  are  listed  in  Table  II.  It  was 
found  that  the  breakdown  shown  as  "Trial  1"  resulted  in  accu¬ 
rate  calculation  of  harmonics  up  to  the  7th.  The  limiting  har¬ 
monic  was  raised  to  the  13th  when  the  size  of  the  azimuthal 
steps  was  decreased  as  indicated  by  "Trial  2"  in  Table  II. 
Finally,  for  hover,  it  was  established  that  harmonics  up  to 
the  20th  will  be  substantially  unaffected  if  the  azimuthal 
increments  are  reduced  to  smaller  values  than  those  associ¬ 
ated  with  "Trial  3".  The  values  shown  as  "calculated"  in  Fig¬ 
ure  7  were  obtained  using  the  "Trial  3"  breakdown. 

In  considering  shortcomings  of  the  theory,  it  must  be  noted 
that  the  prediction  method  has  been  concerned  from  the  outset 
with  only  the  sound  resulting  from  resultant  lift  and  drag 
forces.  Pressure  fields  due  to  blade  thickness,  boundary  lay¬ 
er  turbulance  and  wake  effects  have  all  been  neglected,  and, 
thus,  one  would  expect  underestimation  of  the  actual  sound 
field.  Another  effect,  not  accounted  for  here,  is  suggested 
by  the  fact  that  some  of  the  measured  peaks  in  the  low  har¬ 
monic  range  are  shifted  appreciably  from  integer  multiples  of 
the  fundamental.  Since  the  rotational  noise  "sources"  are 
viewed  as  fixed  on  the  swept  area  in  this  analysis,  such  an 
effect  cannot  be  associated  with  doppler  shifts  due  to  source 
motion.  On  the  other  hand,  much  of  the  sound  is  received  after 
passing  through  the  non-uniform  and  time-varying  vertical  flow 
field  just  beneath  the  rotor.  This  could  cause  frequency 
shifts  and  changes  in  the  radiation  pattern  due  to  refraction. 
It  is  likely,  however,  that  for  radiation  pattern  changes  the 
presence  of  the  ground  reflecting  plane  has  a  more  pronounced 
influence.  A  very  curious  aspect  of  the  experimental  spectrum 
in  Figure  7  is  the  fact  that  the  6th,  9th  and  12th  harmonic 
peaks  do  not  appear  to  have  been  shifted. 

In  spite  of  all  the  uncertainties  in  the  comparison  between 
test  and  calculation  discussed  above,  it  seems  that  for  har¬ 
monics  higher  than  the  3rd,  rotational  noise  is  probably  not 
the  major  source  in  hover.  It  is  worth  noting,  however,  that 
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because  of  the  high  1st  harmonic,  the  calculated  overall  noise  levels 
would  check  reasonably  well  with  measured  overall  levels  even  in 
hover. 

ROTOR  NOISE  IN  FORWARD  FLIGHT 

In  the  machine  computation  of  rotor  noise,  basic  input  data 
were  obtained  from  Ref.  16  and  Ref.  17,  i.e.,  the  data  for  ASL, 
o  ,  3,  8  versus  helicopter  and  flight  velocity.  Table  III 

provides  specific  cross-reference  information.  Tables  V,  VI, 
and  VII  are  exact  print-outs  of  the  airfoil  input  data  used, 
i.e.,  CL,  CD,  CN  versus  a  and  M  for  12  percent,  15  percent  and 

18  percent  symmetric  airfoils,  respectively.  In  view  of  their 
length,  these  three  tables  are  located  immediately  preceding 
the  list  of  References  at  the  end  of  this  report. 

The  effect  of  forward  speed  is  shown  in  Figures  8  through  15 
for  the  H-34  and  HU-1  helicopters.  The  major  influence  of  for¬ 
ward  flight  is  clearly  in  the  higher  harmonics,  as  would  be  ex¬ 
pected.  Since  the  prediction  method  is  limited  to  the  har¬ 
monics  contained  in  the  representation  chosen  for  the  chordwise 
pressure  distributions  shown  in  Figure  3,  it  is  appropriate  to 
examine  the  values  of  the  Fourier  amplitude  coefficients, 

c  ,  d_  and  e_  .  These  are  shown  as  functions  of  harmonic 
m.  mL  m. 

number  for  seven  radial  stations  in  Figures  16  through  21.  The 
"flat"  character  of  these  curves  at  the  more  outboard  stations 
indicates  that  the  pressure  distributions  shown  as  a  function 
of  time  in  Figure  3  are  well  approximated  by  impulses  for  these 
and  higher  values  of  r^/c^* 

The  Fourier  amplitude  coefficients  of  an  impulse  are,  of  course, 
constant  for  all  harmonics  (Ref.  5,  Figure  3).  Deviations 
from  impulse-type  harmonic  content  occur  in  the  inboard  sta¬ 
tions;  it  follows  that  before  eliminating  Eg.  (3)  in  favor  of 
a  constant,  the  relative  importance  of  inboard  versus  outboard 
stations  as  noise  sources  must  be  established.  The  greater 
variations  in  velocity  over  the  inboard  blade  sections  at  high 
forward  speed  would  be  expected  to  emphasize  the  contributions 
of  these  inboard  sections.  Accordingly,  sound  pressure  levels 
were  calculated  for  the  H-34  helicopter,  flying  at  115  knots 
and  a  rotor  speed  of  231  RPM  at  a  field  point  where  Ro**200  feet, 

i|)o=150°,  and  zq— 10.2  feet.  It  had  been  found  that  in  hover, 

the  contribution  of  the  outboard-most  10  percent  of  radius 
could  be  more  than  50  percent  of  the  sound  in  almost  all  har¬ 
monics.  But  at  V*115  knots,  the  important  contributions  come 
from  the  more  inboard  rings.  Accordingly,  revised  integra¬ 
tion  breakdowns,  Trials  4  and  5  in  Table  II,  were  defined. 

The  percentage  contribution  from  the  concentric  rings  of  swept 
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TABLE 
FOR  INPUT 

III 

DATA  FOR 

ASL , 

6,  6 

Helicopter 

Velocity 

Ref. 

ASL 

Paqe 

6/0 

“s 

0 

19 

19 

10 

H-34 

23 

16 

29 

29 

10 

79 

64 

64 

10 

115 

109 

109 

10 

0 

152-153 

317 

5,35,319 

HU-1 

34 

17 

110-111 

269 

5,35,271 

88 

124-125 

285 

5,35,287 

113 

138-139 

301 

5,35,303 
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Ft 0.8 .  MAIN  ROTOR  ROTATIONAL  NOISE 
FOR  H-34  HELICOPTER. 
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Fig. 9  MAIN  ROTOR  ROTATIONAL  NOISE 
FOR  H-34  HELICOPTER. 
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Fig. II.  MAIN  ROTOR  ROTATIONAL  NOISE 
FOR  H-34  HELICOPTER. 
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Flg.12.  MAIN  ROTOR  ROTATIONAL  NOISE 
FOR  HU- 1  HELICOPTER. 
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Fig.13.  MAIN  ROTOR  ROTATIONAL  NOISE 
FOR  HU -I  HELICOPTER. 
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Fig. 14.  MAIN  ROTOR  ROTATIONAL  NOISE 
FOR  HU -I  HELICOPTER. 
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FI«.I5.  MAIN  ROTOR  ROTATIONAL  NOISE 
FOR  HU  — I  HELICOPTER. 
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Fig.  16.  FOURIER  AMPLITUDE  COEFFICIENTS 
FOR  H-34  HELICOPTER. 
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Fig.  17.  FOURIER  AMPLITUDE  COEFFICIENTS 
FOR  HU— I  HELICOPTER. 
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Fig.  19. FOURIER  AMPLITUDE  COEFFICIENTS 
FOR  HU- 1  HELICOPTER. 
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area  are  indicated  in  Table  IV  based  on  a  Trial  4  breakdown. 
Note  that  a  percentage  greater  than  100  means  that  tl:«t  out¬ 
board  rings  had  a  canceling  rather  than  a  reinforcing  effect. 

It  would  be  expected  that  the  higher  the  harmonic  content  of 
the  aerodynamic  section  loading ,  the  more  azimuthal  stations 
might  be  required  to  properly  predict  the  higher  harmonics. 

In  spite  of  this,  no  important  differences  were  noted  between 
Trials  3  and  4,  so  the  Trial  3  breakdown  was  used  for  all  cal¬ 
culations  from  hover  to  maximum  forward  speed. 

It  is  interesting  that  the  increase  in  higher  harmonic  noise 
level  with  forward  speed,  as  shown  in  Figures  8  through  11 
(i.e.,  for  the  H-34  helicopter),  is  rather  gradual.  The  same 
increase  for  the  HU-1  helicopter,  as  shown  in  Figures  12 
through  15,  however,  appears  to  take  place  almost  entirely 
between  hover  and  the  "transitional"  forward  speed  of  34  knots. 
The  cause  of  this  is  most  likely  to  be  in  the  aerodynamic  sec¬ 
tion  loading;  a  comparison  of  this  basic  information  is, 
therefore,  shown  in  Figures  22  and  23  for  the  H-34  and  HU-1, 
respectively.  The  oscillatory  nature  of  these  curves  can  be 
seen  to  increase  markedly  for  the  HU-1  between  V-0  and  V"»34 
knots,  while  there  is  a  gradual  increase  in  the  number  of 
"wiggles"  with  speed  for  the  H-34. 

RADIATION  PATTERNS 

Enough  calculations  have  been  performed  to  allow  sound  radia¬ 
tion  patterns  to  be  examined  for  individual  harmonics.  Such 
a  presentation,  however,  would  entail  a  very  large  number  of 
figures.  Since  field  measurements  are  often  made  using  octave 
band  filters,  it  was  decided  to  present  the  calculated  data 
in  that  form,  in  spite  of  harmonics  falling  on  frequencies 
which  are  borderline  between  one  octave  band  and  another. 
Figures  24  through  27  show  azimuthal  patterns  for  the  H-34  in 
hover  and  at  three  forward  speeds  in  the  first  three  octave 
bands  at  a  radius  of  two  rotor  diameters.  Figures  28  through 
30  show  similar  data  at  a  distance  of  ten  diameters. 

Although  dissymmetries  begin  to  become  appreciable  for  the  H-34 
at  higher  forward  speeds,  there  is  nearly  polar  symmetry.  Only 
at  ten  diameters  and  in  the  first  octave  band  radiation  pattern 
of  the  H-34  at  high  forward  speed,  does  the  pattern  tend  to  be 
lobe-like  on  either  side  of  the  longitudinal  center  line  of  the 
helicopter  and  facing  roughly  aft.  The  existence  of  azimuthal 
asymmetry  apparently  depends  on  both  field  point  radius  and 
speed.  These  changes  with  field  point  radius  are  shown  by 
comparison  of  Figures  25  and  28,  Figures  26  and  29,  and  Figures 
27  and  30,  for  field  points  located  in  a  plane  one  diameter  be¬ 
low  the  helicopter.  The  patterns  might  show  greater  changes 
at  other  values  of  observer  elevation;  these  were  not  investi¬ 
gated. 
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AERODYNAMIC  SECTION  LOADING  Ob/ 
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Fig. 23.  AERODYNAMIC  SECTION  LOADING  ,  HU- 1  HELICOPTER 
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Fig. 26.  RADIATION  PATTERN  IN  OCTAVE  BANDS, 

H—  34  HELICOPTER. 
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Fig.3Q  RADIATION  PATTERN  IN  OCTAVE  BANOS, 

H— 34  HELICOPTER. 

Radiation  patterns  for  the  HU-1  are  available  but,  since  the 
data  indicate  more  symmetry  than  the  patterns  for  the  H-34, 
these  are  not  presented. 

The  results  of  Garrick  and  Watkins5  would  lead  one  to  predict 
that  rotor  radiation  patterns  in  the  zQ  direction  would  show 

greater  variations.  Figures  31  and  32  for  the  H-34  helicopter 
show  that  this  is  indeed  so  for  the  first  three  octaves  at 
Ro/D*2,10  and  4>o*270 •  In  fact,  the  general  character  of  the 

variation  beneath  a  lifting  rotor  is  roughly  the  same  as  that 
of  the  fundamental  behind  a  propeller  in  axisymmetric  flow, 
regardless  of  forward  speed.  Sound  pressure  levels  increase, 
rapidly  in  most  cases,  for  hub  height  downward  and  then  fall 
off  slowly  as  the  zq  distance  is  increased  negatively.  At 

two  radial  diameters  (Figure  31) ,  the  maxima  occur  at  about 
two  diameters  of  negative  elevation;  at  ten  radial  diameters 
(Figure  32)  ,  the  maxima  occur  at  about  five  diameters  of  nega¬ 
tive  elevation. 

Variations  with  radial  distance  are  of  great  interest  for  the 
problem  of  aural  detection.  In  this  instance  the  individual 
frequency  components  are  most  pertinent,  since  attenuation 
due  to  atmospheric  or  terrain  absorption1  is  dependent  on  fre¬ 
quency.  However,  to  reduce  the  data  presented,  the  calculated 
results  are  again  given  by  octave  band.  The  radial  trends  are 
shown  on  Figures  33  and  34  for  the  H-34  at  V«0  and  V-115  knots, 
respectively.  In  both  cases  ^*270°  and  zq/D«-1. 

The  particular  point  of  interest  is  the  distance  beyond  which 
"spherical  spreading"20  can  be  assumed  to  take  place.  Figures 
33  and  34,  therefore,  have,  in  addition  to  the  calculated 
sound  levels,  a  line  which  represents  "spherical  spreading” 
extrapolated  inward  from  a  point  100  diameters  from  the  rotor 
center.  Since  the  correction  for  spherical  spreading  is 

R1 

Adb  -  20  log1Q(^) 

2 

where 

r  «■  distance  from  source  point  where  sound  pressure 
level  is  known 

R2  ■  distance  from  source  to  point  where  one  wishes  to 
know  the  sound  pressure  level, 

and  since  the  abscissa  of  Figures  33  and  34  is  a  log  scale,  the 
curves  of  spherical  spreading  are  straight  lines  with  slopes 
of  -20  decibels  per  decade.  Where  the  calculated  curve  becomes 
tangent  to  and  henceforth  follows  this  line,  then  the  total 
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Fig. 32.  VERTICAL  RADIATION  PATTERN,  H-34 HELICOPTER. 
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Flg.33  CALCULATED  a  SPHERICAL  SPREADING 
FOR  H-34  HELICOPTER. 
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variation  with  radius  can  be  said  to  be  due  to  spherical 
spreading  alone.  If  this  is  the  criterion  for  the  "far  field," 
then  Figure  34  shows  that  the  "far  field"  could  be  one  hundred 
diameters  from  the  source.  Note  in  Figure  34  that  if  a  dis¬ 
tance  of  three  diameters  were  used  as  the  base  for  extrapola¬ 
ting  the  sound  pressure  level  in  the  first  octave  band  to  a 
distance  of  one  hundred  diameters,  as  is  often  suggested,  an 
error  of  10  decibels  could  result. 

TREND  STUDIES 

Compressibility  effects  were  examined  by  comparing  sound  pres¬ 
sure  levels  associated  with  12-percent  and  18- percent  airfoils. 
Since  one  would  expect  larger  compressibility  effects  at  higher 
velocity,  the  trend  is  examined  for  the  H-34  in  the  115-knot 
case.  The  results  are  presented  on  Figure  35.  For  the  field 
point  chosen,  there  is  not  much  change  in  sound  level  between 
the  12-percent  and  18-percent  airfoils. 

The  effect  of  RPM  was  investigated  with  the  H-34  flying  at 
79  knots.  One  notes  that,  in  reality,  there  would  be  a  change 
in  v*V/flR;  this  implies  a  change  in  ASL  variation  with  azi¬ 
muth.  However,  in  order  to  isolate  the  RPM  effect,  such 
changes  in  ASL  are  neglected.  Airfoil  section  characteristics 
are  the  same  for  both  cases;  appropriate  changes  in  0  and  e  are 
made,  however.  The  results  are  given  for  RPM*214  and  245  on 
Figure  36.  For  the  field  point  chosen/  the  only  significant  ad¬ 
verse  effect  of  increasing  RPM  appears  to  be  in  the  fifth  har¬ 
monic. 

The  effect  of  increasing  blade  chord  was  also  investigated, 
again  using  data  from  the  H-34  at  79  knots.  The  comparison  is 
made  of  normal  chord  and  one  1.5  times  normal.  Appropriate 
changes  are  made  in  6  and  6,  assuming  that  blade  weight  and 
pitch  mass  moment  of  inertia  ihcrease  linearly  with  chord. 
Again,  constant  airfoil  section  characteristics  are  assumed. 

The  results  are  presented  on  Figure  37,  where  a  small  but 
fairly  general  decrease  in  sound  level  is  noted  for  the  larger 
chord. 

The  effect  of  changing  the  number  of  blades  from  two  to  three 
was  investigated  for  the  HU-1  in  the  88-knot  case.  The  chord 
and  ASL  are  adjusted  to  preserve  constant  total  lift.  No  ad¬ 
justment  is  required  in  0  and  6,  again  assuming  a  linear  de¬ 
crease  in  blade  weight  with  chord  (which  provides  the  proper 
balance  of  lift  and  centrifugal  moments  to  preserve  0) .  The 
results  are  presented  on  Figure  38;  for  the  field  point  chosen, 
there  is  a  very  significant  decrease  in  sound  level  for  all 
harmonics  due  to  increasing  the  number  of  blades  from  two  to 
three . 
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Fig.  35.  COMPRESSIBILITY  EFFECTS 
FOR  H-34 HELICOPTER. 
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Fig.  36.  EFFE 
FOR  H-341 


R0/D  •  2 
Z0/D  •  -I 
%  •  270* 

V  •  79  KTS 


r-T— l~— T— IT— I  I  I1  I  I  > 

10  15  20 

MIC  NUMBER 


OF  R.P.M. 


RMS  PRESSURE 
(db  ret  to  0.0004  dynes/cm* ) 


2  BLADES  (fund.  frtq.  *10.4  cp») 


I  5  10  15  20 


HARMONIC  NUMBER 


Fig.  38.  EFFECT  OF  NUMBER  OF  BLADES 
FOR  HU -I  HELICOPTER. 
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The  effect  of  impulsive  downwash  was  investigated  by  selec¬ 
tive  change  in  the  aerodynamic  section  loading  of  the  H-34 
in  the  115-knot  case.  The  modification  was  made  over  an  in¬ 
terval  of  azimuth  approximately  20  degrees  on  each  side  of 
100*  and  for  fractional  blade  spans  of  r/R»0. 85 ,0. 90  and 
0.95,  as  shown  graphically  in  Figure  39.  This  modification 
was  suggested  by  the  data  on  Figure  7  in  Ref.  23.  That  is, 
theory  indicates  sharper  variations  than  experiment.  Since 
data  used  in  the  present  work  was  obtained  using  instrumen¬ 
tation  with  a  break-point  of  30  cps,  making  the  ASL  more  like 
the  theory  would  tend  to  correct  for  "rounding  off"  of  sharp 
peaks  due  to  instrumentation  limitations. 

The  results  of  the  ASL  modification  are  shown  on  Figure  40 
which  indicates  that  the  main  effect  of  this  change  is  to  vary 
the  sound  level  more  markedly  from  one  harmonic  to  another,  at 
least  at  the  one  field  point  used  in  the  comparison. 

The  final  trend  studies  of  the  present  investigation  were 
associated  with  setting  8,  and  $  equal  to  zero.  This  last 

change,  of  course,  is  tantamount  to  eliminating  induced  drag. 

The  case  investigated  was  the  H-34  in  hover?  the  field  point 

was  <l>  *270,  R  /D«-l,  z  /D*-l.  The  results  indicate  no  signi- 
o  Q  o 

ficant  effect.  Ref.  10  shows  the  contribution  of  coning  angle 
to  the  sound  field  and  indicates  that  its  maximum  effect  is  in 
the  plane  of  the  disc.  In  the  present  investigation,  with  the 
field  point  one  diameter  below  the  rotor,  the  effect  of  coning 
leaves  the  sound  level  unaffected  to  within  0.5  decibel.  This 
suggests  that  for  typical  helicopter  rotors,  a  greatly  simpli¬ 
fied  program  could  be  used  with  no  significant  loss  of  accuracy. 
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Fig.  40.  EFFECT  OF  IMPULSIVE  DOWNWASH 
FOR  H- 34  HELICOPTER. 
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0*0857 

179*0 

0*0 

-0*09 

0.0170 

0*0903 

179*0 

0*  1 

-0*0900 

0*0170 

0*0903 

179.0 

0*2 

-0*0800 

0.0170 

0*0803 

179.0 

0*3 

-0*0840 

0*0170 

0*0843 

179*0 

0*4 

-0*1000 

0*0170 

0*1003 

179*0 

0*5 

-0*1000 

0*0170 

0*1003 

179*0 

0*6 

-0*1350 

0*01 70 

0*1353 

179.0 

0.7 

-0*1300 

0*0173 

0*1303 

179*0 

1*0 

.-0*0350 

0*1620 

0*0378 

'iTT' 


TABLE  V  -  (Cont'd) 


180*0  0*0 
180*0  0*1 
180*0  0*2 
180*0  0*3 

180*0  0*4 

180*0  0*5 

1 80  *0  0*6 
180*0  0*7 

180*0  1*0 


0*0150  0*0 
0*0150  0*0 
0*0150  0*0 
0*0150  0*0 
0*0150  0*0 
0*0150  .)*0 
0*0450  0*0 
0*0150  0*0 
0*1585  0*0 


o.o 

0.0 

0.0 

•  0066 

0.0 

0*0 

0.1 

0.0 

•  0086 

0.0 

0*0 

0.2 

6.0 

•  0086 

0.0 

0*0 

0.3 

0.0 

•  0086 

0.0 

0*0 

0.4 

0.0 

•  0086 

0.0 

0.0 

0.5 

0.0 

•  0087 

0*0 

0.0 

0.6 

0.0 

•  0087 

0*0 

0.0 

0.7 

0.0 

•  0088 

0.0 

0.0 

i  .0 

0.0 

•  1 068 

0*0 

1.0 

0.0 

.1026 

•  0088 

•  1030 

1  .0 

0.1 

.  1028 

•  0068 

•  1030 

1.0 

0.2 

•  1009 

•  0090 

•  1010 

1  .0 

0.3 

.1006 

•  0090 

•  1006 

1.0 

0.4 

•  1076 

•  0090 

•  1008 

1.0 

0.3 

•  1170 

•  0090 

•  1  171 

1.0 

0.6 

•  1375 

•  0091 

•  1376 

1.0 

0.7 

.  1440 

•  0094 

•  1441 

1.0 

1  .0 

.0250 

•  1036 

•  0266 

2.0 

0.0 

•  2057 

•  0092 

•  2059 

2.0 

0.1 

•  2057 

•  0092 

•  2059 

2.0 

0.2 

*2018 

•  0092 

•  2019 

2.0 

0.3 

•  2012 

.0092 

•  2014 

2.0 

0.4 

•  2152 

•  0092 

•  2154 

2.0 

0.3 

•  2365 

•  0092 

•  2367 

2.0 

0.6 

•  2650 

•  0092 

•  2652 

2.0 

0.7 

•  2790 

•  0103 

•  2792 

2.0 

1.0 

•  0550 

•  1060 

•  0587 

3.0 

0.0 

•  3086 

•  0094 

•  3086 

3.0 

0.1 

•  3086 

••0094 

•  3086 

3.0 

0.2 

•  3026 

•  0094 

•  3027 

3.0 

0.3 

•  3019 

•  0094 

•  3020 

3.0 

0.4 

•  3226 

•  0096 

•  3229 

3.0 

0.5 

•  3510 

•  0096 

•  3510 

3.0 

0.6 

•  3725 

•  0096 

•  3725 

3.0 

0.7 

•  4150 

•  0150 

•  4152 

3.0 

1.0 

•  0925 

•  1074 

•  0980 

4.0 

0.0 

•  4114 

•  0100 

•  411 1 

4.0 

0.1 

•  41  14 

•  0100 

•  411  1 

4.0 

0.2 

•  4036 

•  0100 

•  4033 

4.0 

0.3 

•  4025 

•  0100 

•  4022 

4.0 

0.4 

•  4305 

•  0100 

•  4302 

4.0 

0.5 

•  4605 

•  0100 

•  4601 

4.0 

0.6 

•4630 

•  0108 

•  4626 

4.0 

0.7 

•  5250 

•  0216 

•  5252 

4.0 

1.0 

•  1400 

•  1166 

»  1478 

5.0 

0.0 

•  5143 

•  0106 

•  5133 
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TABLE  VI  -  (Cont'fl 


a 

M 

CL 

CD 

CN 

5*0 

0*  1 

•  5143 

•  0106 

•  5133 

5*0 

0.2 

•  5044 

•  0106 

•  5034 

5.0 

0.3 

.5032 

•  0108 

•  5022 

5*0 

0.4 

.5381 

.0108 

•  5370 

5*0 

0.5 

•  580  0 

•  0108 

•  5787 

5.0 

0.6 

•  5975 

•  0156 

•  5966 

5.0 

0.7 

•  6050 

•  C346 

•  6057 

5*0 

1  .0 

•  2125 

•  1430 

•  2242 

6*0 

0.0 

•  6172 

•  01  18 

•  6150 

6.0 

0.1 

•  6172 

•  01  18 

•  6150 

6*0 

0.2 

•  6053 

•  01  18 

•  6032 

6*0 

0.3 

•  603b 

•  01  18 

•  6017 

6.0 

0.4 

•  6458 

•  01  18 

•  6434 

6*0 

0.5 

•  6895 

•  0127 

•  6863 

6*0 

0.6 

•  6900 

•  0249 

•  6668 

6*0 

0.7 

•  6575 

•  0528 

.6594 

6*0 

1  .0 

•  2650 

•  1710 

•  2814 

7.0 

0.0 

•  7200 

•  0132 

.7162 

7.0 

0.1 

•  7200 

•  0132 

.7162 

7.0 

0.2 

•  7062 

•  0132 

.7025 

7.0 

0.3 

•  7044 

•  0132 

.7008 

7.0 

0.4 

•  7534 

•  0132 

.7494 

7,0 

0.5 

•  7675 

•  0172 

.7639 

7.0 

0.6 

•  7800 

.0356 

.7785 

7.0 

0.7 

•  6810 

•  0750 

•  6851 

7.0 

1  .o 

•  3075 

•  2024 

•  3299 

8.0 

0.0 

•  8100 

•  0146 

•  8042 

8.0 

0.1 

•  8100 

•  0146 

•  8042 

8.0 

0.2 

•  8070 

•  0146 

•  8012 

8.0 

0.3 

•  8050 

•  0146 

•  7992 

8.0 

0.4 

•  8600 

•  0154 

•  8538 

8.0 

0.5 

.8550 

•  C271 

•  8505 

8.0 

0.6 

•  8340 

•  0572 

•  8338 

8.0 

0.7 

•  6975 

•  0905 

.7033 

8.0 

1  .0 

•  3300 

•  2156 

•  3568 

y.o 

0.0 

•  ?100 

•  0162 

•  9013 

9.0 

0.1 

•  9100 

•  0162 

•  9013 

9.0 

0.2 

•  9050 

•  0162 

•  6964 

9.0 

0.3 

•  9050 

•  0162 

•  8964 

9.0 

0*4 

•  9500 

•  0212 

•  9416 

9.0 

0*5 

•  9300 

•  0415 

•  9250 

9.0 

0.6 

•  8550 

•  0772 

•  8566 

9.0 

0.7 

•  7000 

•  1  168 

•  7096 

9.0 

1  *0 

•  3525 

•  2368 

•  3652 

10.0 

0*0 

1  *0125 

•  0182 

1  *0003 

10.0 

0*1 

1*0125 

•  0182 

1  *0003 
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TABLE  VI  -  (Cont'd) 


a 

N 

CL 

CD 

°N 

0*0 

0.2 

1 .0050 

•  0164 

•  9929 

0.0 

0.3 

1 .0000 

•  0195 

•  9882 

o.  u 

0.4 

1 .0200 

.0317 

1*0100 

0.0 

0.5 

.9740 

•  0604 

•  9697 

0*0 

0.6 

•  860b 

•  0985 

•  6645 

0*0 

0.7 

•  7050 

•  1365 

•  7183 

0*0 

1.0 

•  3750 

•  2578 

•  4  141 

1.0 

0.0 

1 .0950 

•  0204 

1*0768 

1  .0 

0.1 

1 .0950 

•  0204 

1 .0788 

1.0 

0.2 

1 • 10U0 

•  0204 

1 .0837 

1  .0 

0.3 

1 .1000 

.0256 

1 .0847 

1  .0 

0.4 

1 • 095U 

.0448 

1 .0634 

1  .0 

0.5 

•  992b 

.0604 

•  9896 

1  .0 

0.6 

•  8450 

•  1 192 

•  6522 

1  .0 

0.7 

•  7030 

•  1568 

•  7204 

1  .0 

1.0 

•  3950 

•  2768 

•  4406 

12.0 

0.0 

1  .2050 

•  0228 

1 *1834 

12.0 

0.1 

1 .2050 

•  0228 

1 • 1834 

L  2.0 

0.2 

1 . 19bO 

•  0232 

1 • 1 694 

12.0 

0.3 

1 • 1 750 

•  0348 

1 • 1 566 

12.0 

0.4 

1.1175 

.0636 

1 • 1063 

12.0 

0.5 

•  9800 

•  1015 

.9797 

12*0 

0.6 

•  8150 

•  1406 

•  6264 

12.0 

0.7 

•  7000 

•  1802 

•  7222 

12.0 

1  .0 

•  4150 

•  2978 

.4679 

13.0 

0.0 

1 .2900 

•  0264 

1 .2629 

13.0 

0.1 

1 .2 900 

•  0264 

1 .2629 

13.0 

0.2 

1 .2675 

•  0301 

1 .2418 

13.0 

0.3 

1 .2050 

•  0496 

1 • 1 853 

13.0 

0.4 

1 .1300 

•  0844 

1.0948 

L  3.0 

0.5 

•  9350 

•  1228 

•  9386 

13.0 

0.6 

•  7775 

•  1616 

•  7939 

13.0 

0.7 

•  6950 

•  2008 

•  7223 

13.0 

1  .0 

•  4375 

•  3178 

•  4978 

14.0 

0.0 

1 .3570 

•  0394 

1 .3262 

14.0 

0.1 

1  .35/0 

•  0394 

1 .3262 

14.0 

0.2 

1  .2970 

•  0472 

1 .2699 

14.0 

0.3 

1 *2025 

•  0750 

1 • 1849 

L4.0 

0.4 

1 .0650 

•  1121 

1 *0605 

14.0 

0.5 

•  8650 

•  1506 

.8952 

14.0 

0.6 

•  7425 

•  1896 

•  7663 

14.0 

0.7 

.6875 

•  2288 

•  7224 

L4.0 

1.0 

.4460 

•  3458 

•  5164 

15.0 

0.0 

1  .3725 

•  0788 

1 .3461 

15.0 

0.1 

1 .3725 

•  0796 

1 *3461 

15.0 

0.2 

1 .2850 

•  0890 

1 .2643 

15.0 

0.3 

1 .1800 

•  1086 

1 • 1679 
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TABLE  VI  -  (Cont'd) 


a  N  CL  CD  °N 


15*0 

0*4 

1*0050 

•  1421 

1  *0075 

15*0 

0*5 

•  8250 

•  181  1 

•  8438 

15.0 

0*6 

•  7100 

•  2214 

•  7431 

15*0 

0.7 

•  6805 

•  2626 

•  7253 

15*0 

1  *0 

•  4815 

•  3846 

•  5646 

165*0 

0*0 

-•591 1 

•  1468 

-•6090 

165*0 

0*  1 

-•591  1 

•  1476 

-•6092 

165*0 

0*2 

-•5498 

•  1574 

-.5718 

165*0 

0*3 

-•5300 

•  1764 

-•5576 

165*0 

0*4 

-•5122 

•  1992 

- *5464 

165*0 

0*5 

-•5131 

.2247 

-•5538 

165*0 

0*6 

-•4636 

•  2540 

-•5329 

165*0 

0*7 

-•4762 

•  2892 

-.5367 

165*0 

1  *0 

-•3978 

•  4154 

-•4917 

166*0 

0*0 

-•6091 

•  1348 

-•6236 

166*0 

0*  1 

-•6091 

•  1350 

-.6237 

166*0 

0*2 

-•571 1 

•  1434 

-•5888 

166*0 

0*3 

-•5488 

•  1606 

-•5714 

166*0 

0*4 

-•5325 

•  1814 

-•5606 

166*0 

0*5 

-.5280 

•  2046 

-•5618 

166*0 

0*6 

-•5016 

•  2301 

-•5426 

166*0 

0.7 

—•4653 

•  2620 

-•5343 

166*0 

1*0 

-•3764 

.3778 

-•4566 

167*0 

0*0 

-•6248 

•  1228 

-•6365 

167*0 

0*  1 

-•6248 

•  1228 

-•6365 

167*0 

0*2 

-•5875 

•  1282 

-•6013 

167*0 

0*3 

-•5626 

•  1438 

-•5806 

167*0 

0*4 

-•5600 

•  1625 

-•5810 

167*0 

0*5 

-.5372 

•  1842 

-•5648 

167*0 

0*6 

-•5132 

•  2088 

-.5471 

167.0 

0*7 

-•4855 

•  2374 

-.5265 

167*0 

1*0 

-•3542 

•  3458 

-•4229 

168*0 

0*0 

-•6226 

•  1  108 

-•6320 

168*0 

0*  1 

-•6226 

•  1  108 

-•6320 

168*0 

0*2 

-•6005 

•  1  1  18 

-•6106 

168*0 

0*3 

-•5785 

•  1240 

-.5916 

168*0 

0*4 

-•5666 

•  1410 

-•5856 

168*0 

0*5 

-•5492 

•  1612 

-.5707 

168*0 

0*6 

-•5232 

•  1854 

-•5503 

168*0 

0*7 

-.4936 

•  2128 

-•5270 

168*0 

1  *0 

-•3372 

•  3198 

-•3963 

169*0 

0*0 

-•6228 

•  0988 

-•6303 

169*0 

0*  1 

-•6228 

•  0988 

-•6303 

169*0 

0*2 

-•6000 

•  0968 

-•6078 

169*0 

0*3 

-•5920 

•  1034 

-•6009 

169*0 

0*4 

-•5722 

•  1 1  76 

-•5841 

169*0 

0*5 

-•5624 

•  1365 

-•5781 
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TABLE  VI  -  (Cont'd) 


a 

N 

CL 

CD 

°N 

169*0 

0*6 

-.5300 

•  1604 

-.5509 

169*0 

0*7 

-•5036 

•  1890 

-•5353 

169*0 

1  *0 

-•3131 

•  2986 

-•3693 

170*0 

0*0 

-•6075 

•  0874 

-.6134 

170*0 

0*  1 

-•6075 

•  0874 

-•6134 

170*0 

0*2 

-.5970 

•  0874 

-•6031 

170*0 

0*3 

-•5900 

•  0880 

-•5963 

170*0 

0*4 

-•5300 

•  0978 

-•5682 

170*0 

0*b 

-.5627 

•  1  141 

-.5740 

170*0 

0*6 

-•5403 

•  1366 

-.5558 

170*0 

0*7 

-.5086 

•  1668 

-•5298 

170*0 

1  *0 

-.3021 

•  2778 

-•3458 

jl  7 1  *0 

0*0 

-•5935 

•  0768 

-•5982 

171.0 

0*  1 

-.5935 

•  0768 

-.5982 

*171*0 

0*2 

-•5668 

•  0766 

-.5915 

171*0 

4 

0*3 

-.5858 

•  0768 

-•5906 

■171*0 

0*4 

-•5842 

•  0808 

-•5896 

171*0 

0*5 

-•5670 

•  0956 

-•5750 

171  *0 

0*6 

-.5472 

•  1160 

-•5586 

;1  7  1  *0 

0.7 

-•5126 

•  1426 

-.5286 

,171*0 

1  *0 

-•2820 

•  2544 

-•3183 

'172*0 

0*0 

-.5532 

•  0658 

-.5569 

172*0 

0*  1 

-•5532 

•  0658 

-.5569 

172*0 

0*2 

-•5700 

•  0658 

-.5736 

172.0 

0*3 

-•5665 

•  0658 

-•5702 

172.0 

0*4 

-.5767 

•  0662 

-•5803 

il  72*0 

i 

0.5 

-.5700 

•  0772 

-.5752 

|l72*0 

0.6 

-.5444 

•  0962 

-•5524 

1172.0 

0.7 

-.51  18 

•  1209 

-•5237 

i  1 72*0 

1.0 

-•2630 

•  2378 

-.2935 

173.0 

0.0 

-•4914 

•  0568 

-•4946 

{173*0 

0.  1 

-•4914 

•  0568 

-.•4946 

173*0 

0.2 

-.5254 

•  0568 

-.5284 

173*0 

0.3 

-•5366 

•  0568 

-•5396' 

173*0 

0*4 

-•5474 

•  0568 

-•5503 

173*0 

0*5 

-.5628 

•  0614 

-•5661 

173*0 

0.6 

-•5400 

•  0767 

-•5453 

173.0 

0.7 

-•5088 

•  0964 

-•5170 

173*0 

1  *0 

-•2409 

•  2188 

-•2658 

174*0 

0*0 

-.4330 

•  0474 

-•4355 

174*0 

0*1 

-.4330 

•  0474 

-•4355 

174*0 

0*2 

-•4659 

•  0474 

-•4683 

174*0 

0*3 

-.4770 

•  0474 

—  <4794 

174*0 

0*4 

-•5223 

•  0474 

-•5244 

174*0 

0*5 

-.5290 

•  0476 

-•531 t 

174*0 

0.6 

-•5378 

•  0578 

-•5409 

174*0 

0.7 

-•4981 

•  0780 

-•5035 
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TABLE  VI-  (Cont'd) 


a 

M 

CL 

CD 

V 

174*0 

1*0 

-•2140 

•  2064 

-•2344 

175*0 

0*0 

-.3697 

•  0388 

-.3717 

175*0 

0*1 

-•3697 

•  0388 

-.3717 

175*0 

0*2 

-•3966 

•  0388 

-•3904 

175*0 

0*3 

-•4025 

•  0380 

-.4044 

1  75*0 

0*4 

-•4520 

•  0383 

-.4537 

175*0 

0*5 

-•4680 

•  0308 

-.4697 

175*0 

0*6 

-•5241 

•  0448 

-•5260 

175*0 

0*7 

-•4840 

•  0635 

-.4877 

175*0 

1  *0 

-•1902 

•  1948 

-•2065 

176*0 

0*0 

-•3114 

•  0308 

-.3127 

176*0 

0*1 

-•3114 

•  0308 

-.3127 

176.0 

0*2 

-.3321 

•  0308 

-.3334 

176*0 

0*3 

-•3280 

•  0300 

-.3293 

176*0 

0*4 

-.3817 

•  0308 

-.3829 

1  76.0 

0*5 

-.3938 

•  0308 

-.3950 

176*0 

0*6 

—•4420 

•  0325 

-.4440 

176*0 

0.7 

-•4568 

•  0478 

-.4590 

176*0 

1*0 

-.1669 

•  1644 

-.1794 

177.0 

0*0 

-.2386 

•  0248 

-.2345 

177.0 

0*1 

-.2336 

•  0248 

-.2345 

177*0 

0.2 

-.2404 

•  0248 

-.2494 

177.0 

0*3 

-.2504 

•  0248 

-.2593 

177*0 

0*4 

-•2913 

•  0248 

-.2922 

177.0 

0*5 

-.3182 

•  0248 

-.3190 

177.0 

0*6 

-.3725 

•  0250 

-.3733 

177.0 

0*7 

-.3724 

•  0367 

-.3739 

177.0 

1*0 

-.1410 

•  1748 

-.1508 

178.0 

'  0*0 

-• 1654 

•  0214 

-• 1660 

170.0 

0*1 

-.1654 

•  0214 

-• 1660 

178.0 

0*2 

-.1506 

•  0214 

-.1660 

170.0 

0*3 

-.1709 

•  0214 

-.1795 

178.0 

0*4 

-.1908 

•  0214 

-.1915 

170.0 

0*5 

-•2150 

•  0214 

-.2164 

178.0 

0*6 

-.2765 

•  0214 

-.2771 

178.0 

0*7 

-•2527 

•  0272 

-.2535 

178*0 

1*0 

-.  1  100 

•  1668 

-• 1 158 

179.0 

0*0 

-•0876 

.0178 

-.0878 

179.0 

0*  1 

-•0876 

•  0178 

-.0678 

179*0 

0*2 

-•0793 

•  0178 

-.0796 

179*0 

0*3 

-•0044 

•  0178 

-.0847 

179*0 

0*4 

-.1004 

•  0178 

-. f007 

179.0 

0.5 

-.1026 

•  0170 

-.1029 

179.0 

0*6 

-• 1435 

•  0178 

-.1430 

179*0 

0.7 

- • 1 356 

•  0212 

-.1360 

J  79*0 

1  *0 

-•0430 

•  1628 

-.0466 
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TABLE  VI  -  (Cont'd) 


a 

M 

CL 

CD 

180.0 

0*0 

-0*0 

•  0)50 

-0*0 

1  BO  . 0 

0*1 

-0*0 

•  0 1 58. 

-0.0 

1  80*0 

0*2 

-0.0 

•  0158 

-0.0 

180.0 

0*3 

-0*0 

•  0158 

-0*0 

180*0 

0*4 

-0*0 

•  0158 

-0*0 

180*0 

0*5 

-0*0 

*0158 

-0*0 

180*0 

0*6 

-0.0 

*0158 

-0*0 

180*0 

0*7 

-0*0 

•  0181 

-0.0 

180*0 

1  *0 

-0.0 

*  1 5f>4 

-0.0 
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TABLE  Vn 

18 -PERCENT  SYMMETRIC  AIRFOIL  INPUT  DATA 


a 

M 

CL 

CD 

0*0 

0.0 

0.0 

•  00*5 

0.0 

0.0 

0.1 

0.0 

•  00*5 

0.0 

u.o 

0.2 

0.0 

•  0095 

0.0 

0*0 

0.3 

0.0 

•  0095 

0.0 

0*0 

0.4 

0.0 

.00*5 

0.0 

0*0 

0.5 

u.o 

•  00*5 

0.0 

0*0 

0.6 

0.0 

•  00*5 

0.0 

0*0 

0.7 

0.0 

•  0096 

0.0 

0*0 

1 .0 

0.0 

•  1075 

0.0 

1  .0 

0.0 

•  1000 

.00*7 

•  1002 

1  •  o 

0.1 

•  1000 

•  0097 

.1002 

1  •  0 

0.2 

•  1000 

•  00*8 

•  1002 

1  *0 

0.3 

•  1000 

•  00*8 

•  1002 

1*0 

0.4 

•  1081 

.00*8 

•  1083 

1  *0 

O.b 

•  1200 

.00*9 

•  1202 

1*0 

0.6 

•  1600 

.00*9 

•  1601 

1.0 

0.7 

.1500 

•  0105 

•  1502 

1  .0 

1.0 

•  0300 

•  1  1  10 

•  031* 

2.0 

0.0 

•  2000 

•  0100 

•  2002 

2.0 

0.1 

.2000 

•  0100 

.2002 

2.0 

0.2 

•  2000 

•  0100 

•  2002 

2.0 

0.3 

•  2000 

•  0100 

•  2002 

2.0 

0.4 

•  2162 

•  0100 

•  2164 

2.0 

0.5 

•  2450 

•  0101 

•  2452 

2.0 

0.6 

•  3000 

•  0101 

•  3002 

2.0 

0.7 

.2620 

•  0120 

•  2822 

2.  J 

1  .0 

•  0700 

•  1140 

•  073* 

3.0 

0.0 

.3000 

.0103 

•  300 1 

3.0 

0.  1 

•  3000 

.0103 

•  3001 

3.0 

0.2 

.3000 

•  0103 

•  3001 

3.0 

0.3 

•  3000 

•  0103 

•  3001 

3.0 

0.4 

•  3243 

•  0104 

•  3244 

3.0 

0.5 

.3600 

*0104 

•  3601 

3.0 

0.6 

.4000 

.0104 

•  4000 

3.0 

0.7 

.4400 

*0183 

•  4404 

3.0 

1.0 

•  1100 

.1165 

•  1  15* 

4.0 

0.0 

.4000 

•  0109 

•  3**6 

4.0 

0.1 

.4000 

.0109 

•  3996 

4.0 

0.2 

.4000 

•  0109 

•  3998 

4.0 

0.3 

.4000 

•  0109 

•  3998 

4.0 

C  .4 

.4324 

•  0109 

•  4321 

4.0 

0.5 

•  4650 

•  0109 

•  4646 

4.0 

0.6 

.4660 

*0123 

•  4657 

4.0 

0.7 

•  5500 

.0253 

•  5504 

4.0 

1.0 

•  1500 

•  1223 

•  1582 

5.0 

0.0 

•5000 

•  01.15 

.4991 
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TABLE  VII  -  (Cont'd) 


a 

M 

CL 

CD 

°N 

5  •  0 

0.1 

•  5000 

•  01  15 

•  4991 

5*0 

0.2 

•  5000 

•  0115 

•  4991 

5.U 

0.3 

•  5000 

.01  16 

•  4991 

5*  O' 

0.4 

•  5405 

•  01  16 

•  5395 

5.U 

0.5 

•  5900 

.0116 

•  5688 

5*0 

0.6 

•  6250 

•  0184 

•  6242 

5.0 

0.7 

•  6100 

•  0372 

•  6109 

5.0 

1  .0 

•  2150 

.1440 

•  2267 

6.0 

0.0 

•  6000 

•  0127 

•  5960 

6.0 

0.1 

•  6000 

•  0127 

.5980 

6.0 

0.2 

•  6000 

•  0127 

•  5980 

6.0 

0.3 

•  6000 

•  0127 

•  5960 

6.0 

0.4 

•  6480 

•  0127 

•  6464 

6.0 

0.5 

•  7000 

•  0144 

•  6977 

6.0 

0.6 

•  7000 

•  0284 

•  6991 

6.0 

0.7 

•  6550 

•  0538 

•  6570 

6.0 

1  .0 

•  2700 

•  1718 

•  2865 

7*0 

0.0 

•  7000 

•  0140 

•  6965 

7.0 

0.1 

•  7000 

•  0140 

•  6965 

7.0 

0.2 

•  7000 

•  0140 

•  6965 

7.0 

0.3 

•  7000 

•  0141 

•  6965 

7.0 

0.4 

•  7567 

•  0141 

•  7528 

7.0 

0.5 

•  7850 

•  0199 

•  7816 

7.0 

0.6 

•  7800 

.0398 

•  7790 

7.0 

0.7 

•  6820 

•  0760 

•  6662 

7.0 

1  .0 

•  3150 

•  2032 

•  3374 

6.0 

0.0 

•  8000 

•  0155 

•  7944 

3.0 

0.1 

•  8000 

•  0155 

•  7944 

8.0 

0.2 

•  8000 

•  0155 

•  7944 

tt.O 

0.3 

•  8000 

•  0155 

•  7944 

ti.O 

0.4 

•  8700 

•  0165 

•  8636 

6.0 

0.5 

•  8700 

•  0322 

•  8660 

8.0 

0.6 

•  8330 

•  0600 

.8332 

8.0 

0.7 

•  7000 

•  0665 

•  7052 

8.0 

1  .0 

•  3400 

•  2164 

•  3666 

9.0 

0.0 

•  9000 

•  0170 

•  6916 

9.0 

0.1 

•  9000 

•  0170 

•  6916 

9.0 

0.2 

•  9000 

•  0171 

•  6916 

9.0 

0.3 

•  9000 

•  0171 

•  8916 

9.0 

0.4 

•  9600 

•  0240 

•  9519 

9.0 

0.5 

•  9300 

•  0460 

•  9257 

9.0 

0.6 

•  8550 

•  0784 

•  8567 

9.0 

0.7 

•  7050 

•  1  175 

•  7147 

9.0 

1  .0 

•  3650 

•  2377 

•  3977 

1  0.0 

0.0 

1*0000 

•  0191 

•  9661 

10.0 

0.1 

1 *0000 

•  0191 

•  9861 
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TABLE  VII  -  (Cont'd) 


0 

M 

CL 

CD 

°N 

10*0 

0.2 

1  .OOOO 

.0192 

•  9881 

10.0 

0.3 

1 .0000 

•  0206 

•  9684 

1  0*0 

0.4 

1 #0200 

•  0378 

1  *01  1  1 

10.0 

0.5 

.9700 

•  0635 

•  9663 

10*0 

0.6 

•  8610 

•  0995 

•  865 2 

1  o.O 

0.7 

•  7100 

•  1395 

•  7234 

10.0 

1  *0 

.3900 

•  2587 

•  4290 

1  1  .0 

0.0 

1 . 1 ooo 

•  02  1 2 

1 • 0638 

1  1  .0 

0.1 

1  . 1000 

•  0212 

1 *0836 

1  1  .0 

0*2 

1 . 10OO 

•  0212 

1 .0838 

11.0 

0.3 

1  • 1000 

•  0265 

1 *0852 

1 1  .u 

0.4 

1 .0800 

•  0505 

1 *0698 

1  1  .0 

0.5 

•  9950 

•  0818 

•99233 

1  1  .0 

0.6 

•  8450 

•  1200 

•  6524 

1  1  .0 

0.7 

.  7080 

•  1596 

•  7254 

1  1*0 

1  *0 

•  4100 

•  2777 

•  4555 

12.0 

0.0 

1  *2100 

•  0237 

1 • 1 685 

12.0 

0.1 

1  .2100 

•  0237 

1 • 1885 

12*0 

0.2 

1  *2000 

•  0245 

1 • 1 789 

12.0 

0.3 

1*1720 

.0395 

1 • 1546 

12.0 

0.4 

1*1 150 

•  0667 

1 • 1045 

12*0 

0.5 

•  9750 

•  1030 

.9751 

12.0 

0.6 

•  8200 

•  1420 

•  8610 

12.0 

0.7 

•  7050 

•  1814 

•  7273 

12*0 

1  .0 

•  4300 

•2987 

•  4832 

13*0 

0.0 

W3000 

•  0272 

1 *2728 

13.0 

0.1 

1  .3000 

•  0272 

1 *2728 

13*0 

0.2 

1  *2750 

•  0332 

1 .2496 

13*0 

0.3 

1.2G0J 

•  0549 

1 • 1816 

13*0 

0.4 

1  *1 1O0 

•  0858 

1 • 1009 

13*0 

0.5 

•  930 0 

•  1240 

•  9341 

13*0 

0.6 

•  7900 

•  1626 

•  8063 

13*0 

0.7 

•  7000 

•  2015 

•  7274 

13*0 

liO 

•  4550 

•  3187 

•  5150 

14.0 

0*0 

1 *3550 

•  0402 

1 .3245 

14.0 

0.1 

#1  .3550 

•  0402 

1 .3245 

14.0 

0*2 

*1  *2950 

.•0532 

1 *2694 

14.0 

0.3 

W20O0 

* *0785 

1 • 1634 

14.0 

0*4 

1  *0500 

•  1132 

1 *0462 

14*0 

0*5 

•  8900 

•  1520 

•  9003 

14.0 

0*6 

•  7600 

•  1906 

•  7835 

14.0 

0.7 

•  6950 

•  2295 

•  7299 

14*0 

wo 

•  4800 

•  3467 

•  5496 

lb.0 

0*0 

1  *3750 

•0797 

1 *3468 

15*0 

0*1 

W3750 

•  0808 

1 *3491 

15*0 

0*2 

1  *2800 

•  0950 

1 *2610 

15*0 

0*3 

1*1800 

•  1  157 

1 • 1697 
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TABLE  VII  -  (Cont'd) 


a 

M 

CL 

CD 

=H 

15.0 

0.4 

.9800 

•  1440 

.9839 

15.0 

0.5 

.8300 

•  1820 

•  8488 

15. U 

0.6 

•  7300 

•  2224 

•  7627 

15.0 

0.7 

•  6850 

•  2632 

•  7298 

lb.O 

1  .O 

•  5030 

•  3854 

.5856 

165.0 

0.0 

—.5922 

.1477 

-.6103 

165.0 

0.1 

-.5922 

•  1493 

-.6107 

165.0 

0.2 

-.54  n 

•  1652 

-.5718 

165.0 

0.3 

-.5300 

•  1858 

-.5600 

165.0 

0.4 

-.4995 

•  2085 

-5364 

165.0 

0.5 

-.5162 

•  2326 

-.5588 

165.0 

0.6 

-.4972 

.2594 

-.5474 

165.0 

0.7 

-.4813 

•  2918 

-.5404 

165.0 

1  *o 

-.4155 

.4162 

-.5091 

166.0 

0.0 

-.6082 

•  1357 

-.6230 

166.0 

0.1 

-.6082 

•  1360 

-.6230 

166.0 

0.2 

-.5702 

•  1486 

-.5892 

166.0 

0.3 

-.5477 

•  1679 

-.5721 

166.0 

0.4 

-.5250 

•  1900 

-.5554 

166.0 

0.5 

-.5310 

.2125 

-.5666 

166.0 

0.6 

-.5137 

•  2352 

-.5553 

166.0 

0.7 

-.4906 

•  2644 

-.5400 

166.0 

1  .0 

-.3927 

•  3787 

-.4727 

167.0 

0.0 

-.6297 

•  1237 

-.6414 

167.0 

0.1 

-.6297 

•  1237 

—.64 1 4 

167.0 

0.2 

-.5910 

•  1310 

-.6053 

167.0 

0.3 

-.5603 

•  1482 

-.5793 

167.0 

0.4 

-.5600 

•  1677 

-.5834 

167.0 

0.5 

-.5343 

.1895 

-.5632 

167.0 

0.6 

-.5215 

.2120 

-.5558 

167.0 

0.7 

-.4890 

•  2392 

-.5303 

167.0 

1  «o 

-.3683 

•  3467 

-.4369 

168.0 

0.0 

-•6252 

.1117 

-.6348 

168.0 

0.1 

-.6252 

•  1  1  17 

-.6348 

168.0 

0.2 

-.6030 

•  1  130 

-.6133 

168.0 

0.3 

-.5770 

.  1280 

-.5910 

168.0 

0.4 

-.5675 

•  1468 

-.5856 

168.0 

0.5 

-.5464 

•  1675 

-.5693 

168.0 

0.6 

-.5264 

.  1900 

-.5544 

168.0 

0.7 

-.4971 

•  2156 

-.531 1 

168.0 

1  .0 

-.3494 

•  3207 

-.4084 

169.0 

0.0 

-.6257 

•  0997 

-.6*332 

169.0 

0.1 

-.6257 

•  0997 

-.6332 

169.0 

0.2 

-.6000 

•  0997 

-.6080 

169.0 

0.3 

-.5920 

•  1068 

-.6015 

169.0 

0.4 

-.5643 

•  1230 

-.5774 

169.0 

0.5 

-.5638 

•  1435 

-.5808 
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TABLE  VII  -  (Cont'd) 


169*0 

169.0 

169*0 

170*0 

170*0 

170*0 

170*0 

170.0 

170.0 

170.0 

170*0 

170*0 

171*0 

171*0 

171.0 

171*0 

171.0 

171.0 

171.0 

171.0 

171*0 

172*0 

172*0 

172.0 

172.0 

172*0 

172.0 

172.0 

172.0 

172.0 

173.0 

173.0 

173.0 

173.0 

173.0 

173.0 

173.0 

173.0 

173.0 

174.0 

174.0 

174.0 

174.0 

174.0 

174.0 

174.0 

174.0 


0.6 
0.7 
1.0 
0.0 
0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
1.0 
0.0 
0.  1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
1.0 
0.0 
O.  1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
1.0 
0.0 
0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
1.0 
0.0 
0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 


-.5300 
-.5122 
-.3302 
-.6000 
—  •  "jOOO 
-•£941 
-.5900 
-.5600 
-.5604 
-.5406 
-.5122 
-.3142 
-.5870 
-.5870 
-.5835 
-.5625 
-.5903 
-.5670 
-.5440 
-.5151 
-.2920 
-.5463 
-.5463 
-.5650 
-.5630 
-.5634 
-.5800 
-.5437 
-.5137 
•—•2709 
—  .477  T 
-.477/ 
-.5206 
-.5333 
-.5499 
-.5757 
-.5400 
-.5095 
-.2466 
-.4209 
-.4209 
—.4018 
-.4741 
-.5246 
-.5371 
-.5456 
-.4962 


.1675 

•  1922 

•  2997 

•  0882 
•  0882 
•  0882 

•  0895 

•  1020 

•  1205 

•  1426 

•  1695 
.2787 
.0777 

•  0777 
.0777 
.0777 

•  0834 

•  1000 

•  1204 

•  1453 

•  2552 
.0667 

•  0667 

•  0667 

•  0667 
.0673 

•  0808 

•  1005 

•  1244 
.2387 

•  0577 

•  0577 

•  0577 

•  0577 

•  0577 

•  0643 

•  0814 

•  1033 

•  2197 

•  0482 

•  0482 

•  0482 

•  0482 

•  0482 

•  0488 

•  0636 
*0657 


-.5522 

-.5395 

-.3813 

-.6062 

-.6062 

-.6004 

-.5966 

-.5889 

-.5728 

-.5572 

-.5339 

-.3578 

-.5919 

-.5919 

-.5865 

-.5874 

-.5961 

-.5757 

-.5561 

-.5315 

-.3253 

-.5503 

-.5503 

-.5668 

-.5668 

-.5871 

-.5856 

-.5524 

-.5260 

-.3015 

-.4812 

-.4812 

-.5240 

-.5364 

-.5528 

-.5792 

-.5459 

-.5163 

-.2717 

-.4236 

-.4236 

-.4643 

-.4765 

-.5266 

-.5393 

-.5493 

-.5024 


93 


TABLE  VIX  -  (Cont'd) 


a 

M 

CL 

CD 

=K 

174*0 

1  *o 

“•2161 

•  2072 

-•2386 

175*0 

0*0 

-•3594 

•  0397 

-•3615 

175*0 

0*1 

-.3594 

•  0397 

-•3615 

175.0 

0*2 

-•3931 

•  0397 

-•3951 

175*0 

0*3 

-•4000 

•  0397 

-•4019 

175*0 

0*4 

-•4540 

.0397 

-•4557 

175*0 

0*5 

-.4761 

•  0397 

-•4777 

175*0 

0*6 

-.5482 

•  0486 

-•5503 

175*0 

0*7 

-•4880 

•  0705 

-•4923 

175*0 

1  *0 

-. 1925 

•  1957 

-•2088 

176*0 

0*0 

-•302  t 

•  0317 

-•3042 

176*0 

0*1 

-.302  7 

.0317 

-•3042 

176*0 

0*2 

-.3292 

•  0317 

-•3306 

176*0 

0*3 

-•3259 

•  0317 

-.3273 

176*0 

0*4 

-.3634 

•  0317 

-.3847 

1  76*0 

0*5 

-.397  / 

•  0317 

-•3969 

176*0 

0*6 

-•4457 

•  0340 

-.4470 

176*0 

0*7 

-.4785 

•  0532 

-•4810 

176.0 

1  *o 

-. 1788 

•  1852 

-• 1913 

177.0 

0*0 

-.2271 

•  0257 

-•2281 

177*0 

0.1 

-•2271 

•  0257 

-•2281 

177.0 

0*2 

-•2469 

.0257 

-•2479 

177.0 

0*3 

-•2567 

•  0257 

-.2577 

177.0 

0*4 

-•2926 

•  0257 

-•2935 

177.0 

0.5 

-•3263 

•  0257 

-•3272 

177.0 

0*6 

-•4000 

•  0260 

-•4008 

177*0 
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